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GENERAL INTRODUCTION 


It has long been recognized by light aircraft manufacturers that the 
technical consideration which exerts the largest influence on sales is the 
performance of the aircraft. If an aircraft can climb higher, fly further, 
cruise faster, and land and take off on shorter runways than other aircraft 
for the same payload and price, then it will usually sell better than its 
competition. A substantial portion of the engineering effort expended on 
a new design is therefore devoted to estimating the performance improvement 
resulting from a given change in configuration or powerplant. What one would 
like to be able to do is to suggest those changes which result in optimum 
performance for the price. 

The National Aeronautics and Space Administration undertook the present 
study to assist the light aircraft industry and aeronautical education in 
general. It was felt that the task of selecting the best performing con- 
figuration for the price would be greatly facilitated if the pertinent 
research results of the last 35 years were readily at hand in easi ly-usable 
forms. A computer program employing a collection and arrangement of those 
methods and data most applicable to the estimation of light aircraft per- 
formance would appear to satisfy these requirements. The present work seeks 
to provide these programs along with a detailed review of the methods used 
and some worked-out examples. The work is thus analogous to previous studies 
(Refs. I, 2) which approached the prediction of riding and handling qualities 
of light aircraft in a similar fashion. 

As will become evident from the literature review following, performance 
estimation can be treated at three levels of sophistication. The first, which 
may be termed static or point performance, is concerned with the maximum 
values for certain parameters such as speed, rate of climb, etc. assuming 
that nothing changes with time. The equations expressing the power required 
for level flight as a function of speed and altitude are algebraic and 
therefore fairly easy to evaluate. The variation in maximum power available 
as a function of flight speed and altitude can be evaluated point by point, 
if not analytically. The region between these two data sets on the speed- 
power plane is that for which steady flight is possible. Thus, by finding 
intersections of the two functions one has maximum and minimum speed while 
the maximum difference between the two curves is a measure of the maximum 
rate of climb. Other preformance parameters are calculated with similar 
directness. The calculations are usually performed graphically because an 
analytical solution requires the extraction of the roots of a fourth or 
higher order polynomial, a laborious procedure if done by hand. It will be 
recognized also that determining the effect of a change in configuration or 
power plant involves many calculations if one wishes to see the effect at all 
weights and operating altitudes. 

Until about 20 years ago all aircraft manufacturers used these perform- 
ance estimation techniques which were first developed in the early 1930's. 

The methods are generally as reliable as the quality of the input lift, drag, 
and thrust data. The results are easily interpreted and can be checked in 



flight through appropriate tests. With the advent of the modern digital 
computer, the major airframe manufacturers began to give consideration to 
more sophisticated means of describing the manner in which an airplane 
performs. They recognized for example that the time required to reach a 
given altitude could be minimized by varying the speed as the altitude 
increases and that the range could be increased on some aircraft by allowing 
the altitude to increase as the fuel is burned. In other words, the path 
over which the aircraft flies determines the performance of the vehicle. 

Hence they began to integrate the differential equations which describe 
aircraft motion with various types of control inputs to see what paths are 
produced. The digital computer permits one to investigate a large number of 
cases quickly and relatively inexpensively. Through a trial and error 
process one can get a good indication of how to fly a particular mission to 
obtai'n optimum results. 

This level of sophistication is in common use among the large airframe 
constructors today. It has long been recognized, however, that an optimum 
path obtained in this fashion cannot be shown to be an optimum in the mathe- 
matical sense. The development of mathematically optimum flight paths has 
been a subject of theoretical research for at least a hundred years. Solutions 
of several simple problems have been obtained but a general procedure that is 
successful in a large number of cases has thus far eluded formulation. 

The purpose of this report is to critically review methods available for 
estimating most aspects of light aircraft performance at all three levels of 
sophistication and to render those methods which are regarded as most accu- 
rate into fast, easy-to-use forms employing a digital computer. Through this 
device it is hoped that light aircraft designers can investigate a wider 
range of parameters economically in their search for improved performance in 
their vehicles. The programs and explanations are written at such a level 
that they should be readily intelligible to recent B.S. graduates. 

As written here, the vehicle lift, drag, and thrust terms in the per- 
formance equations are represented by implicit functions. To obtain numer- 
ical solutions, explicit functions are required. These the program obtains 
by making rather general fits of user-supplied data. Unfortunately, it was 
net possible within the scope of the present work to eliminate the require- 
ment for the user to supply these data. It would have been desirable to ask 
the user to specify only the aircraft geometry and the power plant and pro- 
peller characteristics and to have the program compute the lift, drag, and 
thrust characteristics needed for the performance computation. 

The work begins with a review of the pertinent literature of the past 40 
years. Estimation techniques based on the point performance concept are then 
developed. These techniques have been programmed for computer solution. The 
use of this program is then explained and some sample results for a typical 
light aircraft are given. 

The next section treats the path performance concept. Again, an easily- 
used computer program has been developed to perform the computations. Its use 
and basis are explained and typical results are provided. 
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An appendix provides a detailed derivation of the equations for path 
performance while listings of the Fortran IV programs used to compute point 
and path performance are given in two add itiona I appendices. 

Other appendices present programs for fitting power curves, lift-drag 
curves, the basis for the integration technique used on the path performance 
equations, and a more detailed discussion of the nature of the fuel-flow- 
power relationship in piston engines. 

The reader will perhaps note the absence of a reference to the standard 
text by Perkins and Hage (John Wiley 1949) and the failure to follow the 
nomenclature of this text which is by now fairly standard. However, it 
seemed that because the equations selected for computer solution are really 
simplifications of those used in stability analysis, the notation should 
follow that common in stability analysis. Some modifications in this 
view were found to be necessary in order to accomodate the more general 
drag polar used in the present work. It is hoped that these departures 
from common usage will not prove too disconcerting. 
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LITERATURE REVIEW 


Point Performance 


"General Formulas and Charts for the Calculation of Airplane Performance", 
TR-408, by Oswald (Ref. 3) and "General Airplane Performance", TR-654, by 
Rockefeller (Ref. 4), published in 1932 and 1939 respectively, represent the 
state of the art in the prediction of point performance. Oswald’s work 
presents a series of performance charts for airplanes equipped with modern 
unsupercharged engines and fixed-pitch metal propellers; these charts yield 
the performance characteristics (maximum level flight speed, maximum rate 
of climb, service ceiling, absolute ceiling, etc.) as a function of the 
parasite drag loading, effective span loading, and thrust horsepower loading. 
Oswald later extended his analysis to include the case of supercharged 
engines (Ref. 5) while White and Martin (Ref. 6) made a similar analysis for 
the case of constant-speed propellers with no supercharging. In each of 
the analyses mentioned above special assumptions were made regarding the 
variation of engine power with altitude and engine speed and the variation 
of propulsive efficiency with altitude and air speed. These assumptions 
along with the assumption of a parabolic drag polar are necessary to obtain 
a problem which is tractable by hand solution techniques or in closed form. 

Rockefeller decided that with new engine and propeller developments it 
would be desirable to attack the problem in a more general manner in order 
to obtain a method of performance calculation basically independent of the 
particular engine-propeller combination but readily adaptive to any type. 

Thus, he developed the equations for the analysis of the performance of an 
ideal airplane — an airplane for which the thrust power is independent of 
speed, the parasite drag is constant, and the lift coefficient has an 
infinite maximum value — in order that the charts developed for use in prac- 
tical calculations would for the most part apply to any type of engine- 
propeller combination and system of control, the only additional material 
required consisting of the actual engine and propel ler curves for the 
propulsion unit. Rockefeller also presented his results graphically as 
performance charts. 

Accurate prediction of point performance characteristics requires 
reliable information on the power the aircraft can put into the airstream. 

For propeller-driven aircraft NACA TR-640 (Ref. 7) and WR L-286 (Ref. 8) 
present propeller data obtained from aerodynamic wind tunnel tests. The 
data is presented as a series of four design charts for each propeller 
tested; these charts have been the standard NASA format since 1929 (see 
Appendix F, Figure (F-1) for an example). Although its basic intent was 
to reveal the effects of changes in solidity resulting either from increasing 
the number of blades or from increasing the blade width, TR-640 is probably 
more widely known for its outline of the procedures required to compute 
the propel ler thrust from the propel ler design charts. A step-by-step 
procedure for calculating the power available is given in Appendix F of 
the present work along with a set of propeller design charts for the 
R.A.F. 6 two blade propeller. 
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The value of knowing the static performance characteristics is voiced 
by Thompson in Reference 9. He mentioned that one of the most perplexing 
guessing games in cross country flying is choosing the most favorable 
altitude and true airspeed for cruising flight. As a means of solving the 
cruising dilemma for level flight with a light airplane, normally operating 
engine, and constant speed propeller he suggested that: 

(1) The high speed dash should be made at near sea level at maximum 
power. 

(2) Normal cruising at 65-75$ power should be made at the highest 
altitude at which these powers are available using full throttle 
and normal cruising RPM. 

(3) Maximum range airspeed should be 1.4 to 2.0 times the flaps up 
stall speed depending on aerodynamic cleanness. 

(4) Range is independent of altitude if airspeed is maintained at 
correct best range speed for each altitude. 

(5) For best range at higher airspeeds, the optimum altitude is 
progressively higher. 

(6) In moderate headwinds, the speed for maximum range should be 
increased about 10$. 

(7) For maximum endurance, the airplane should be flown between 20 
and 30 percent above flaps up stall speed, depending upon where 
minimum power is required to sustain’ level flight. 

The suggestions given by Thompson are generally in good agreement with 
the results obtained from a point performance analysis of the Cessna 182 
(see the section on Examples of Point Performance Calculation). Similar 
agreement was also found using a path performance analysis when flying 
near the angle of attack for best lift to drag ratio. These analyses were 
made using the point and path performance programs presented in Appendices 
C and D respectively. 

In recent years new interest has arisen in improving the performance 
of light aircraft. As noted in Reference 10 the basic technology and 
configurations of most of the present light airplane fleet were developed 
before the advent of the high speed computer, jet transport, high lift 
technology, advanced stability and control analysis methods, analytical 
descriptions of handling qualities, and greatly improved wind tunnel testing 
techniques. Since this advanced technology has not been widely applied to 
light aircraft, they have not kept pace with the improvements achieved by 
commercial airliners. Roskam and Kohl man found by parametric variation that 
aerodynamic design modifications can be made to improve significantly the 
performance of light aircraft. They used a relatively simple computer 
program to evaluate the speed for best range, maximum level flight speed, 
specific range, maximum rate of climb, and speed for maximum rate of climb 
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in terms of the predicted lift and drag coefficients resulting from specific 
geometric modifications. 

Accurate prediction of static performance requires good estimates of 
the lift and drag coefficients as a starting point. An ideal procedure 
for obtaining suitable values of these coefficients would require that 
one specify only the body coordinates, speed, and altitude of the airplane 
to obtain in a precise fashion both the lift and drag coefficients as 
functions of angle of attack; unfortunately, such a procedure is not as 
yet available. Historically, the drag coefficient has been much more 
difficult to estimate accurately than the lift coefficient. Reference 11 
is an example of a sophisticated method for obtaining aerodynamic charac- 
teristics of multi -component a i rfoi I s— ai rfoi Is with leading or trailing 
edge high lift devices — in subsonic viscous flows. The calculated aerody- 
namic characteristics include pressure distribution, lift, pitching-moment, 
and skin friction drag up to incipent separation on any component. The 
characteristics are obtained from a computer program written for either the 
UN I VAC 1108 or the CDC 6600 computer which requires the inputs of freestream 
conditions and the airfoil geometry. Similar techniques are needed to 
handle the complete wing-body-tail combination. 

Two recent works should be helpful in predicting the drag coefficients 
of light aircraft. Roskam in Reference 12 presents two methods for computing 
drag polars of airplanes at subsonic Mach numbers. The first method models 
the drag polar by Cq = CDo + C^/neAR and then sums the zero-lift drag 
coefficients (usually from wind tunnel data) of each individual component 
of the aircraft. For a more detailed and accurate drag prediction Roskam 
suggests a second method. This method employs formulae and charts to estimate 
the zero- lift drag coefficient of the wing-body combination, the horizontal 
tail, and the vertical tail as functions of thickness to chord ratio. The 
drag of the wing-body due to lift Is considered to be a function of wing 
drag due to lift and body drag due to angle of attack; a procedure is also 
given for estimating incremental drag coefficient due to miscellaneous 
components such as windshields, nacelles, flaps, etc. 

Reference 13 by Wolowicz and Yancey which describes methods for estimating 
the longitudinal aerodynamic characteristics of light, twin-engine, propeller- 
driven airplanes, presents a method for estimating the drag coefficient 
very similar to the second method given by Roskam and discussed above. Also 
presented are methods for obtaining lift coefficients of the wing, fuselage, 
horizontal and vertical tails, and interference effects. Most of the 
methods mentioned above require the use of charts or the manual evaluation 
of formulas to obtain the lift and drag. A computer program to speed up 
and mechanize this process would materially simplify point and path per- 
formance estimation. 

Discussions of point performance are incomplete without some considera- 
tion of take-off and landing. NACA TR-450 by Walter S. Diehl (Ref. 14) is 
concerned with the development of a method suitable for routine take-off 
calculations which is reasonably simple without neglecting any important 
variables (See the section on Take-Off and Landing Performance for the 


7 



general equation). While the method presented in the Technical Report is 
intended as a practical approximation to a difficult problem, Diehl believed 
that a more accurate method probably would have no significance in view of 
the crude state of the lift, drag, and thrust data. Diehl therefore reduced 
the ground run formula to S = K s v|/(T^/W) where V s is the take-off speed, 

Tj is the initial net acceleration force, W is the take-off weight, and K s 
is a coefficient depending only on the ratio of initial to final net 
acceleration force. A relation to estimate the time required to take-off 
is also given. 

A more exact approach to take-off and landing performance is given in 
Reference 15 prepared by Boeing Aircraft Company. The Boeing report gives 
a derivation of the basic take-off and landing equation leaving it in 
integral form. Provided the thrust, lift, drag, and load factor due to 
rotation during the approach are known, a numerical integration technique 
can thus be used to evaluate the take-off ground run, time to lift-off, 
ground distance while climbing to 50 feet, time to climb to 50 feet, 
ground distance from 50 feet to touchdown, and ground run after touchdown. 

A detailed discussion of both Reference 14 and Reference 15 is included in 
the Take-Off and Landing Performance section. 


Path Performance 

Aircraft can often exceed the equilibrium value of maximum speed, 
altitude, rate of climb, etc. during periods of accelerated flight. The 
designer seeking the ultimate in vehicle performance will wish to devise 
trajectories which maximize particular parameters of interest. One method 
for doing this, which is growing in popularity with the capability to 
manipulate and apply it, is the technique of specifying schedules of two 
control parameters and determining the motion resulting therefrom. The 
expression "growing in popularity", however, should be used somewhat 
advisedly. One sees indications of the use of such techniques in the 
literature and private conversations with industry people also point in 
the same direction, but specific solutions or calculation procedures are 
noticeably absent. The results given in Reference 16 represent elementary 
forms of such procedures. 

It is quite probably that several computer programs are currently in 
use which will compute the trajectory of an aircraft by integrating the 
first order ordinary differential equations of motion (the integration 
becomes possibly only when some of the unknown parameters such as power, 
lift and drag, velocity, etc. are specified as functions of time so as to 
yield the same number of unknowns as equations). Apparently, the programs 
are either classified or used only for in-house work by the companies who 
developed them because none have been described in the open literature. 

They are thus unavailable to the academic community or the light aircraft 
industry. For example Reference 17 indicates the existence of a landing 
analysis digital computer program developed by the Air Force Flight Dynamics 
Laboratory. This program evolved from a need for comprehensive, quantitative 
analysis of aircraft take-off and landing character i st ics . Similar programs 
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no doubt exist at most of the major aircraft manufacturing companies. Although 
the light aircraft designer may not need an extremely sophisticated procedure 
for performance prediction, he should be provided with a procedure which 
permits him to realize some of the performance improvements resulting from 
modern technology and which frees him from complete reliance on the basic 
design charts of the 1930’s. 


Optimum Performance Paths 


Mathematicians have long been concerned with finding the trajectory which 
optimizes a particular performance criterion. One may cite the classical 
brachistochrone problem — find the shape of a wire along which a frictionless 
bead will move under the influence of its own weight from the origin to some 
other point in minimum time — as a simple example. In attempting to apply 
variational techniques to the flight of powered aircraft, however, one finds 
that the additional degrees of freedom present in a realistic mathematical 
model lead to an almost intractable problem. At first sight, the determina- 
tion of an optimal range trajectory could appear to be capable of treatment 
as a classical Mayer problem (see Reference 18 for statement) but mathemati- 
cal difficulties, apparently encountered by all who have attempted this 
approach, have proven insurmountable. Those who have employed basically 
this approach with success have considered more restricted problems such as 
unpowered flight (Ref. 19). 

Within the last 10 years the subject has received intense study because 
of its applicability to the trajectories of spacecraft. However, one prin- 
cipal feature of many of these analyses — the absence of aerodynamic drag — 
makes them inapplicable to aircraft use. In general, while the aircraft 
problems treated by the newer methods of dynamic programming and the 
Pontryagin Maximum Principle have improved in realism, the techniques are 
still too complex and too restricted for general computational use. The 
reader interested in the details of the newer, more successful mathematical 
methods is directed to References 20, 21, 22, 23, and 24. The last is a 
particularly good treatment of the subject. 

It may be pointed out that while true optimum trajectories can be 
obtained only with the methods cited, practically speaking there are many 
near-optimum trajectories which differ little from the optimum in terms of 
the value of the performance criterion. Some of these near-optimum 
trajectories can usually be found without excessive difficulty by interative 
use of the path performance techniques discussed previously. 


9 



POINT PERFORMANCE 


INTRODUCTION 


The process of predicting an aircraft’s static or point performance 
reduces to an investigation of the condition of the flight path assuming that 
the dependent variables (except for h) of the performance equations do not 
change with time. If for the general performance equations (Appendix B) the 
derivatives of the dependent variables are neglected while assuming n to be 
a new dependent variable, a set of non-linear algebraic equations are obtained. 
The values of the dependent variables required to obtain an optimum flight 
condition (i.e. maximum rate of climb, maximum level flight speed, etc.) can 
be found by applying the Ordinary Theory of Maxima and Minima. 

The work of Oswald (Ref. 3), Rockefeller (Ref. 4), and White and 
Martin (Ref. 6) as noted earlier represents the state of the art in static 
performance prediction. Each of these works choose a parabolic relationship 
between lift and drag in order to obtain a tractable problem while Reference 3 
employs in addition some special assumptions regarding the variation of power 
with velocity for different types of propellers. It is felt, however, that 
a significant improvement can be made in static performance prediction (1) by 
using a drag polar which is more general than the conventional parabolic 
polar and (2) by permitting the user to specify only several points on a 
curve of maximum power available versus velocity rather than the functional 
form of the power-velocity relationship. This course has been followed In 
the present work. 

The static performance equations given herein were developed from those 
in Appendix B by applying the Theory of Maxima and Minima. In addition, a 
general power versus velocity curve and a drag polar of the form 

C D = k 1 + k 2 C L + k 3 C L 4 

have been employed. Note that the parabolic polar is a special case of the 
above polar with and k 4 equal to zero. The quantities which the analysis 
considers as known are: 

C D (C L ) = kj + k 2 C* + k^C l 4 where the user specifies the k’s. Alternately 

the k's may be determined from experimental 
data using the procedure in Appendix E. 

S = wing area on which Cq and C|_ are based, 

W = ai rp lane weight, 

h = altitude at which the optimum characteristics are desired, 

P(V) = points on the maximum power available versus velocity curve 
at specified reference altitude. 

Velocity is the unknown and one desires the velocity for which a particular 
flight characteristic is optimum. An additional restriction imposed on the 
general equations of motion (Appendix B) for the present analysis is that 
cos y 2 1.0. Should the reader be interested in flight conditions with y 
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greater than twelve to fifteen degrees he is referred to the program in 
Appendix D which integrates the general equations of motion without making 
the assumption of small values of flight path angle. Once the data defined 
above is provided the following quantities can be calculated by use of the 
digital computer program listed in Appendix C: 

(1) maximum and minimum level flight speed at an altitude, 

(2) speed for maximum climb angle and maximum climb angle at an altitude, 

(3) speed for minimum power (maximum endurance) and minimum power at 
an altitude, 

(4) classical speed for maximum range, 

(5) service ceiling and the velocity at service ceiling, 

(6) absolute ceiling and the velocity at absolute ceiling, 

(7) maximum rate of climb schedule from altitude hj to altitude , 

(8) most economical rate of climb schedule from altitude hj to 
altitude 

(9) maximum rate of climb, power available and power required versus 
velocity at an altitude. 

Note that because of the generalized nature of the drag and power relations 
employed here, obtaining these optimal quantities usually requires either the 
solution of a pseudo-polynomial equation having non-integer exponents with 
velocity as the unknown, or the solution of a pair of the equations with 
both velocity and altitude as unknowns. 

Discussed in detail in succeeding sections are: 

(1) the derivation of the optimal static performance equations, 

(2) a description of the computerization procedure required to find 
the solution of these equations, 

(3) examples of the static performance calculations for two light 
ai rcraft, 

(4) a discussion of techniques used to calculate landing and take-off 
performance.* 

* This discussion is based on two methods contained in the literature and 
may therefore differ from other methods currently being used. A general 
computer program to evaluate take-off and landing performance has been 
omitted because of the difficulty encountered in estimating correct values 
of C|_, Cq, velocity, and power in these two flight modes. 

12 


■ 


DERIVATION OF THE POINT PERFORMANCE EQUATIONS 


The point or static performance equations are derived by requiring that 
the acceleration terms be zero. With this restriction the equations developed 
in Appendix B become 


x = (V cos y)At + x Q 

(D 

fi = V sin y 

(2) 

aP Sp o c d v2 ° 

WV ' 9 T" W " 9 Sln Y = 0 

(3) 

Sp C.V 2 o 

g — - — g cos y = 0 

(4) 


2 W 

where 

a = f (h) = ( 1 - 6.86 x lO^h) 4,26 . 


(Note that the x equation is given in its integrated form.) In the above set 
of equations h will be treated as a variable separate from by this device 
the equation becomes a system of four algebraic equations in nine unknowns 
(x, V, y, h, h, P, W, C|_, Cq). It is customary to specify aircraft weight, 
power as a function of altitude and velocity, C|_ and Cp as a function of 
angle of attack, and one other parameter in order to make the system solvable. 


Now, if one were to seek information on the flight path angle possible 
at a given altitude he would first express Cq as a general function of C^, 

C D = C D 0 + Xrf + kjC^, 

and then from equation (4) write Cj_ as 

r = 2W cos Y 
L Sp Q 0V 2 ' 

and finally substitute these expressions into equation (3) to yield: 


(5) 


(6) 


P _ Sp. 

W ” 


[c + K ( 2W cos k ( 2W COS Y )M- V sin y = 0 . 
2W [ D o 1 Sp G aV 2 2 Sp Q aV 2 J 1 


(7) 


Note that if k 2 = 0 and kj = 1/(iTeAR) the drag polar in equation (5) reduces 
to the familiar parabolic form. Equation (7) expresses y in terms of V when 
P, h, and W are given. To find the maximum y it is necessary to find the 
value of velocity for which dy/dV = 0. The reader will recognize that this 
is not easily done in closed form. Because of this difficulty, the physics of 
the situation are invoked to reduce the mathematical complexity. Most general 
aviation aircraft do not have sufficient power to climb at angles in excess 
of 12° or 13°; thus, the assumption that cos y = 1.0 is never in error by 
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more than two or three percent. Further, it is difficult to determine Cq q 
and e with errors of less than four or five percent. It may therefore 
be argued that no additional error of significance is introduced by taking 
cos y = 1.0. With this assumption equation (7) becomes 


P 

W 


2W 


C n + k i (' ■ " W --> 2 + k ( — — — ) k 3 1 - V sin y 
Do 1 Cn rr V 2 2 e~ rr \, 2 1 


Sp 0 aV' 


Sp G aV' 


3 - 


= 0. 


( 8 ) 


Finally, through the use of equation (2) In can be written as 

\c 


. P Sp o aV3 


W 


2W 


C 


+ kl ( - 2W - ) 2 + k 2 ( — W -) k 3 

L-'o ' n _ _\/2 ^ c ^ ^\/ 2 


Sp 0 aV' 


Sp 0 aV‘ 




(9) 


This is the fundamental point performance equation which expresses the rate 
of climb as a function of speed when P C V , h ) , W, and h are given. 


Maximum and Minimum Level Flight Speed 


The maximum and minimum level flight speeds are found by setting h = 0 
in the fundamental equation and solving for V in 


p _ Sp 
W 


•CV 3 r 

2 C D 

2W |_ 


+ k , ( — — . ) 2 + k ? ( 
° 1 sp 0 av 2 


2W_) k 3 


2 'sp 0 aV 2 


> 


( 10) 


Maximum Rate of Climb 


The maximum rate of climb occurs at that speed for which 


dh = 1 dP Zld 
dV W dV 2W 


[■ 

o gv3 r 

2W L 


Cn. + k, (^Ur) 2 + k 2 (^U) k 3 


D 0 T 1 
4k 


Sp o aV< 


S P 0 aV ' 


5 


1 ( -2W_>. . 


2k 3 k 2 


SP 0 a 




(2k 3 + 1 ) Sp Q a 


5 


( 11 ) 


is zero. This speed, substituted into equation (9) gives the maximum rate 
of climb. 


Maximum Climb Angle 

An equation for the climb angle can be obtained by noticing that the 
flight path angle has already been assumed to be small. Hence, fi=VsinY~ 
V y> °r Y = h/V. Division of equation (9) by V therefore yields the desired 
relation. The speed for the maximum value of y is then found by setting 
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2W 


dx = 0 = -dE.J E_ 

dV dV WV WV 2 


S Po qV 

w 



k 1 


2W 


S Po aV2 


•) 2 + k 2 ( 


SpoOV' 


’] 


Sp aV 2 r 4k! 9W 

^ 2W ^ 2 

2W V s SP Q a 


2k 3 k 2 

v (2k 3 + 1) 



and solving for V. This speed is then substituted into the equation for y 
to obtain the maximum climb angle. 


Service Ceiling and Absolute Ceiling 


The service ceiling is that value of h for which the maximum value of 
h has decreased # to 100 feet per minute. The absolute ceiling is that value 
of h for which h = 0 at the speed for maximum rate of climb. One must solve 
two equations simultaneously in order to find the value of V for which a is 
a minimum when h = 100 ft/min or h = 0. The two equations are repeated below 
as equations (13) and (14). 


h - r: + 


P Sp^aV 3 


s 


W 2W ^ l/ q 
2 


II- / 2W \ 2 ii- , 2W M 
+ k 1 «w'2> + k 9 ( c ~ ZUT ) = 0 


'1 v Sp 0 a\F' 2 Sp 0 aV 


1 dP 3S Po CTV ^ 


W d V 


( 


Crs + ( 


2W )* + k. (, 


2« \ T "1 " 2 


2W r) kj 


]■ 

] 


sp 0 av 3 r 

2W [" 


1_ k ( _2W_ )2 _ 2k 3 k 2 

V 5 1 S P 0 P (2k 3 + 1) ' S Po° 


2W , k 3 i _ 


(~=-> ^ = 0 
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h is then found from the expression 


(1 - 


h = 


1 

.4.26 


■)10 e 


6.86 


(feet) . 


(13) 


(14) 


Maximum Range Speed 


Fuel consumption in propeller driven aircraft is generally directly 
proportional to the power required. Thus for maximum range at y = 0, one 
should fly at the speed for which the ratio of power required per unit speed 
is a minimum 


d(P/V) = 
dV 


0 = Sp Q aV 


+ 


f£s 2 !l 

2 



k 1 


(^W_,2 + k ( ^W_) k 3 
Sp G aV 2 Sp Q aV 2 


Sp O 0 


... 2k 3 k 2 (_2W_) k 3 

y (2k 3 + 1 ) Sp oa 


( 15 ) 


15 


and solving for V gives the velocity for maximum range. For jet aircraft, 
the fuel consumption is more appropriately taken to be directly related to 
the thrust required. 


Maximum Endurance 


For maximum endurance In propeller-driven aircraft one is interested in 
flying at the speed for which the power required is a minimum or, the speed 
which satisfies the equation: 


dy = 0 = 2 s Po aV 

s p 0 ° v3 r 



2W 


‘Sp 0 aV 


t)‘ 


k 2 ( 


2W 

Sp Q aV 


2 


V 


2k 3 k 2 
(2k 3 + 1) 


A 

Sp Q o 



(16) 


The minimum power can then be found by substituting this value of velocity 
into equation (9) with li = 0. 


Mini mum Time to Climb 


The minimum time to climb is the shortest time required to climb from 
one altitude to another altitude. It can be expressed In intergral form by 



(17) 


where the velocity for maximum rate of climb is found from equation (11) and 
the maximum rate of climb is then evaluated using this velocity in equation 
(9). 


Most Economical Climb 


The most economical climb is that climb technique which will move an 
aircraft from h^ to h2 while using the least fuel, df. Since df = - dW and 
W = - cP for propeller aircraft, the following procedure may be used to 
minimize df/dh. 


df dW dW/dt 

dh = “ dh dh/dt 


cP. 


W 

fi-FTW 

w w 


WcP 

r~-~Dv • 


The above expression will have its minimum value when V is a solution to 


dV 


= 0 . 
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Thus 


9 


Hf dfx WcP , 

d dh J d P - DV 


WcP 

dP 

dD 1 

dV ~ 0 “ dV 

" P - DV ~ 

CM 

> 

Q 

1 

Q_ 

dV 

- V dV - D J 


or. 


Since 


— DV - PV — - PD = 0 . 
d V dV 


D = | p n V 2 S a 


% 


4W z k 


2W 


(Sp 0 ) 2 a 2 V 4 k 2 ( sp 


-) k 3 (— — ) k 3 
aV 2 
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then 


dD 4W 2 kj « 

~ = p Q SaCn V 2- -U- 

dV Ko S p 0 a V 3 


•f 


(k 3 - 1)(p 0 Sak 2 )( 


_2W 

SP ( 


f’1 


v (2k 3 - 


( 18 ) 

(19) 

( 20 ) 


Substituting (19) and (20) into (18) yields an equation which can be solved 
for the velocity for most economical climb. The power used to solve equation 
(18) is the maximum power available to the aircraft; this is evident from 
the following sketch of df/dh versus P (here only power greater than DV is 
considered since this is the minimum power required for climbing flight). 


Wc P . df 
P - DV d h 



Power 


Figure 1. Minimum value of df/dh as a function of power. 


df/dh has its minimum value when the power is maximum. Thus, one should fly 
at the maximum possible power and the velocity given by equation (18) for 
most econom i ca I climb. 
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COMPUTERIZATION PROCEDURE FOR POINT PERFORMANCE EQUATIONS 


A computer program* to evaluate the static performance of an aircraft 
was written in Fortran IV for use on the IBM 370-165 computer. The procedure 
employed can best be described by considering the three major portions into 
which the programming task was divided: 

(1) the expression of the maximum power available in a general func- 
tional form having a smooth, continuous first derivative given a 
set of experimental or calculated maximum power-velocity values, 


( 2 ) 


the application of a least-squares-distance curve-fitting technique 
to fit lift and drag data with a general drag polar of the form 


C D = k, + ♦ k 3 C L k 4 


(3) the utilization of the method of regula falsi (false position) to 
find the roots of the pseudo-polynomials in velocity derived in 
the previous section. 


These three divisions will now be described in more detail. 


Maximum Power Available 


For propeller driven aircraft the maximum power available is a function 
of both velocity and altitude. However, if a maximum power available versus 
velocity curve is known at some reference altitude the power available at 
some other altitude may be obtained from the reference curve by means of a 
multiplicative correction factor (Ref. 3) depending solely on altitude. 
Denote the power available at velocity V and altitude h by P gv (V,h) and the 
power available at velocity V and the reference altitude h re f by P re f(V). 
Then, for an unsupercharged engine 


P av (v » h) 


Prefix 


a - 0.165 } 

a ref * 0*165 


vhere: 


a = p/p Q = ratio of the density of air at altitude h to 
“ the density at sea level, 


0 re f = p re f/p 0 = ratio of the density of air at reference 
altitude to the density at sea level. 


The variation of the density ratio may be approximated by a(h) = (1.0 - 
6.86 x 10 -6 h) 4,26 where h is in feet. 


This program is found in Appendix C along with instruction for its use. 
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For a supercharged engine the variation of power available with altitude 
is small up to some fixed altitude which depends on the particular super- 
charger. Theref6re, in this work the altitude variation of power available 
for a supercharged engine is assumed to be zero, i.e. P av (V,h) = P re f(V), 
and the reference altitude is taken to be sea level. 

The general nature of the power available versus velocity curve at the 
reference altitude considering both fixed pitch and constant speed propeller 
aircraft makes it very difficult to fix this curve with a constant coeffi- 
cient polynomial having a smooth, we I l-behaved first derivative over the 
entire velocity range. Since the derivative of the power available appears 
in the static performance equations, a smooth first derivative is a necessity 
for obtaining valid solutions to these equations. The spline technique has 
been shown to give the best mathematical fit for a set of data points, and 
the cubic spline is the simplest fit which provides a well-behaved first 
derivative. For this reason the maximum power available data points were 
fitted using the cubic spline given in Reference 25 with modified end 
conditions. 

The power available at some velocity V at the reference altitude was 
denoted above by P re f(V). Let the data points from the reference power 
available versus velocity curve be denoted by (P re f)i» Vj, i = 1,2,...,N, 
where N is the number of data points. Then the cubic spline fit is of the 
form 

P ref (V) = Cj V 3 + C 2 V 2 + C 3 V t C 4 

J J j j 

where: 

Vj £ V V J+1 , and j = 1,2,...,N-1 


Note that for each of the N-1 intervals the cubic's coefficients depend on 
the interval (Vj, Vj + j) in which the velocity V lies, but they are constant 
in a given interval. It is this variation of the coefficients that gives 
the spline fit its remarkable curve fitting properties. 

To develop a suitable maximum-power-avai lable-at-any-altitude relation- 
ship for use in the static performance calculations, the computer program* 
(Appendix C) requires that the user supply only a set of power available 
versus velocity data points at a reference altitude, the reference altitude, 
and a control parameter denoting whether or not the engine is supercharged. 

A technique for estimating these data points is given in Appendix F. 


Lift-Drag Curve Fitting Technique 

Most previous estimation techniques for static performance have relied 
on a conventional parabolic drag polar. Because some performance parameters 


* The program in Appendix D also used this procedure to obtain maximum 
power avai lab le. 
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are evaluated at relatively high angles of attack where the drag is 
frequently greater than predicted by a parabolic fit to high speed data, the 
use of the parabol ic polar for such computations leads to erroneous results. 
Accordingly, the static performance analysis in the preceeding section 
permits the use of a general drag polar of the form 


C 


D 


+ k 2 + k 3 



Note that this general drag polar includes the parabolic polar as a special 
case, i.e. k^ = k^ = 0. 

If lift and drag data, preferably up to CL max , are available, the coef- 
ficients of the general drag polar can be obtained with the program given in 
Appendix E. This program, a modification of the one given in Reference 26, 
uses a I east-squa res-d i stance technique to fit the data, i.e. it minimizes 
the sum of the squares of the perpendicular distances fr6m the data point 
to the fitted curve. This type of least squares technique is desirable 
because the drag coefficient versus lift coefficient curve has regions of 
both small and large slopes. 


The program gives the user the option of the following four particular 
forms of this general drag polar: 


(1) C D = k, + ko C? + k, C,' 4 




s 2 


'3 U L 


(2) C = C + k„ C? + k, C, 4 


D 


D 


2 L 


3 L 


( 3 ) C D = k 1 + k 3 C L 4 

(4) C D = c Dq + k ; c*« 


(k 2 = 0) 


(k 2 = 0) 


In cases (1) and (3) all coefficients are varied in the fitting process, and 
ki may not be the actual zero-lift drag coefficient. In cases (2) and (4) 
the user specifies Cq 0 , the zero- lift drag coefficient, and it is not varied 
in the fitting process.* 


Solution of the Pseudo-Polynomials 

The determination of each performance parameter requires the solution 
of a pseudo-polynomial** in velocity, except the cases of service ceiling 
and absolute ceiling where the simultaneous solution of two coupled 


* Cases (3) and (4) are the drag polar forms used in the program which 
integrates the equations of motion (Appendix D). 

** The word pseudo-polynomial refers to a polynomial which has the unknown 
raised to powers which are not necessarily integers. 
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pseudo-polynomials, one in velocity with its coefficients depending on 
altitude and one in altitude with its coefficients depending on velocity, 
are required. 

For those cases where the solution of a single pseudo-polynomial is 
required, the method of regula falsi (false position) is used (Ref. 27). 

The velocity range of the power available versus velocity curve is searched 
by increments until a sign change of the pseudo-polynomial occurs. The 
false position method is then used to obtain the zero of the pseudo- 
polynomial to within a specified tolerance. Note that since the spline 
curve fit (see above section) of the power available versus velocity gives 
a cubic polynomial with different coefficients from each velocity interval, 
the coefficients of the various pseudo-polynomials will vary as the entire 
velocity range is searched. 

In the case of service or absolute ceilings the method of false position 
is used in conjunction with an iteration between the velocity polynomial and 
the altitude polynomial. This iterative procedure is as follows: 

(1) A service or absolute ceiling is assumed, and the velocity pseudo- 
polynomial is solved by vegula falsi. 

(2) With this velocity root the altitude pseudo-polynomial is solved 
by regula falsi for a new value of the service or absolute 

cei I ing. 

(3) With this new altitude step (1) is repeated. Iteration continues 
between step (1) and step (2) until convergence on both velocity 
and altitude is achieved to within a specified tolerance. 
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EXAMPLES OF POINT PERFORMANCE CALCULATIONS 


In order to evaluate the static performance program, two single engine 
light aircraft were investigated; both the description and the results of 
these two test cases are given below. 

The performance of the Cessna 182 (Figure 2) was evaluated using a 
conventional parabolic drag polar and a polar obtained by curve fitting lift 
and drag coefficients obtained from Cessna Aircraft through personal commu- 
nication. The reference wing area for the lift and drag coefficients was 
174 square feet while the basic weight was 2650 pounds. The maximum power 
available curve was obtained using the procedure outlined in Appendix F. 

The engine, a Continental Model 0-470-R, had a power rating of 230 BHP at 
2600 RPM. The maximum engine speed for continuous operation was assumed 
to be 2400 RPM. The propeller had a R.A.F. 6 section with a diameter of 
seven feet. 

Tables 1 and 2 may be used to compare the performance of the Cessna 182 
with the two different drag polars; these tables present the performance 
characteristics for the parabolic and the fitted drag polars respectively. 

The major difference is seen when comparing the minimum level flight speeds. 
The analysis with a parabolic polar gives a much lower value of minimum level 
flight speed than does the fitted polar; this is caused by the error en- 
countered when using a parabolic polar at large lift coefficients (high 
angles of attack). 

The program may be easily adapted to find variations of static perfor- 
mance parameters for different values of aircraft weight and altitude. The 
power available and required curves at sea level, 8000 feet, and 16000 feet 
are presented in Figure 3 for the Cessna 182 with the fitted drag polar. 

The variations in maximum and minimum level flight speeds and maximum rate 
of climb with weight and altitude are given in Figure 4. Figure 5 has also 
been included to indicate how the changes in weight affect the service and 
absolute cei I ings. 

Because of the availability of lift and drag data (Ref. 28) the per- 
formance of the Navi on was evaluated using the static performance program. 
Several points from the wind tunnel data were used to find a general drag 
polar*, while a 285 BHP engine rate at 2900 RPM was used at the power plant. 
The results of this analysis are given in Table 3. 


* The sample output given in Appendix E contains the actual data for the 
Navi on aircraft. 
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Figure 2. Three view of the Cessna 182 







FONER AVAILABLE VS. VELOCITY 
REFERENCE ALTITUDE - 0.0 FEET 

FA IFT-LB S/SEC I VIFT/SECI 

0.0 0.0 

0.291500 05 0.273)00 02 

0.524700 05 0.54*700 02 

0.A9960D 05 0.020000 02 

0.81*200 05 0. I 09330 03 

0.874500 05 0.13**70 03 

0.909480 05 0.1*4000 03 

0.9444*0 05 0.191330 03 

0.958*50 05 0.218*70 03 

0.9*7780 05 0.24*000 03 

0.973*10 05 0.273330 03 

0.979440 05 0.300*70 03 

0.994700 05 0.328000 03 

0.994700 05 0.355330 03 

0.994700 05 0.382**0 03 


AIRCRAFT CHARACTERISTICS 

CO • 0.2*9000-01 • 0.44024D-01*Cl**2 ♦ 0.0 «CL*« 0.0 

MINS AREA * 0 . I 74 COO 03 SQ.FT HEIGHT - 0.2*5000 04 LBS 


STATIC PERFORMANCE AT AN ALTITUDE • 0.0 FT 

N1TH MINIMUM TIME ANO MOST ECONOMICAL CLIMB SCHEDULES TQ A FINAL ALTITUDE - 0.100000 OS FT 


MINIMUM LEVEL FLIGHT SPEED - 0.383470 02 FT /SEC 
LIFT COEFFICIENT - 0.8703*0 01 DRAG COEFFICIENT • 0.33*180 01 

MAXIMUM LEVEL FLIGHT SPEED « 0.253850 03 FT /SEC 
LIFT COEFFICIENT - 0.198610 00 ORAG COEFFICIENT « 0.28*370-01 


MAXIMUM CLIMB ANGLE * 0.128580 02 OEG 
VELOCITY FOR MAXIMUM CLIMB ANGLE * 0.853090 02 FT/SEC 
LIFT COEFFICIENT • C. 175860 01 DRAG COEFFICIENT • 0.1*3050 00 


VaOCITY FOR MAXIMUM ENDURANCE * 0.972250 02 FT/SEC 
POWER FOR MAXIMUM ENDURANCE • 0.2047*0 05 FT-LBS/SEC 
LIFT COEFFICIENT - 0.135390 01 ORAG COEFFICIENT - 0.107*00 00 


VELOCITY FOR CLASSICAL MAXIMUM RANGE • 0.1279*0 03 FT/SEC 
LIFT COEFFICIENT « 0.781660 00 ORAG COEFFICIENT • 0.530000-01 


SERVICE CEILING » 0.228310 05 FT 
VELOCITY AT SERVICE CEILING - 0.150450 03 FT/SEC 
LIFT COEFFICIENT - 0.11*810 01 ORAG COEFFICIENT - 0.8*9700-01 

ABSOLUTE CElLtNG - 0.2487*0 05 FT 
VELOCITY AT ABSOLUTE CEILING - 0.154700 03 FT/SEC 
LIFT COEFFICIENT - 0.118670 01 ORAG COEFFICIENT - 0.8889*0-01 


MAXIMUM RATE OF CLIMB SCHEDULE FROM 0.0 FT TO 0.100000 05 FT 


HIFTI 


R/CIFT/SEC) 

VIFT/SECI 

PIFT-LBS/SECI 

CL 


CO 

TISECI 


0.0 


0.23**40 

02 

0.127*50 

03 

0.859940 

05 

0.78S38D 

00 

0.540550-01 

0.0 


0.500000 

03 

0.230980 

02 

0.1278*0 

03 

0.845340 

05 

0. 794)90 

00 

0.54*820-01 

0.213880 

02 

O.IOOOOD 

04 

0.225370 

02 

0.128070 

03 

0.030900 

05 

0.8034*0 

00 

0.55)190-01 

0.433040 

02 

0.150000 

04 

0.219800 

02 

0.1283Q0 

03 

0.01**30 

05 

0.812570 

00 

0.559*80-01 

0.657710 

02 

0.200000 

04 

0.214270 

02 

0.128520 

03 

0.602520 

05 

0.821740 

00 

0.5**270-01 

0.888130 

02 

0.250000 

04 

0.208790 

02 

0.1287*0 

03 

0.788580 

05 

0.830950 

00 

0.572970-01 

0.112450 

03 

0.300000 

04 

0.203350 

02 

0.129010 

03 

0.774790 

05 

0.840200 

00 

0.579780-01 

0.13*720 

03 

0.350000 

04 

0.197950 

02 

0.129260 

03 

0.7*11*0 

05 

0.849490 

00 

0.58**90-01 

0.1*1*50 

03 

0.400000 

04 

0.192*00 

02 

0.129530 

03 

0.7*7*90 

05 

0.858810 

00 

0.593700-01 

0.187250 

03 

0.450000 

04 

0.187290 

02 

0.129800 

03 

0.734370 

05 

0.8*81*0 

00 

0.600810-01 

0.21)580 

03 

0.500000 

04 

0.182020 

02 

0.130090 

03 

0.721200 

05 

0.877530 

00 

0.608010-01 

0.240*10 

03 

0.550000 

04 

0.17*790 

02 

0.130380 

03 

0.708190 

05 

0.88*920 

00 

0.615300-01 

0.2*8540 

03 

0.600000 

04 

0.171*00 

02 

0.130*80 

03 

0. *95330 

05 

0.89*320 

00 

0.622*80-01 

0.297250 

03 

0. *50000 

04 

0.1**450 

02 

0.131000 

03 

0.682*20 

05 

0.905730 

00 

0.630 150-01 

0.32*840 

03 

0.700000 

04 

0.1*1340 

02 

0.131330 

03 

0.6700*0 

05 

0.915140 00 

0.637*90-01 

0.357350 

03 

0.750000 

04 

0.15*270 

02 

0.131*70 

03 

0.657*40 

05 

0.924540 

00 

0.645310-01 

0.388850 

03 

0.800000 

04 

0.151230 

02 

0.132020 

0) 

0. *45370 

05 

0.9)39)0 

00 

0.652990-01 

0.421)80 

03 

0.850000 

04 

0.146230 

02 

0.132380 

03 

0.63)240 

05 

0.943300 

00 

0.6*0730-01 

0.455000 

03 

0.900000 

04 

0.141270 

02 

0.1327*0 

03 

0.6212*0 

05 

0.952*50 

00 

0.6*8530-01 

0.489800 

03 

0.950000 

04 

0.13*350 

02 

0.133150 

03 

0. *09420 

05 

0.9*19*0 

00 

0.67*380-01 

0.5258)0 

03 

O.IOOOOD 

05 

0.1314*0 

02 

0.13)550 

03 

0.597710 

05 

0.971220 

00 

0.684270-01 

0.4*3180 

03 


Table 1 . Performance of Cessna 182 with parabolic drag polar. 
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POWER AVAILABLE VS* VELOCITY 
REFERENCE ALTITUDE - 0.0 FEET 

PA<FT-LBS/SEC) VIFT/SECI 

0.0 0.0 

0.291500 05 0.273500 02 

0.324700 05 0. 596700 02 

0.699600 05 0.820000 02 

0.816200 05 0.109330 03 

0.874500 05 0.136670 03 

0.909480 05 0.164 000 0 3 

0.944460 05 0.191330 03 

0.958650 05 0.216670 03 

0.967780 05 0.246000 03 

0.973610 05 0.273330 03 

0.979440 05 0.300670 03 

0.994700 05 0.328000 03 

0.994700 05 0.355330 03 

0.994700 05 0.382660 03 


AIRCRAFT CHARACTER! STICS 

CO * 0.268800-01 • 0.56242D-OIKl**2 • 0. 17751D-0l*CL** 0.65000D 01 
VINO AREA - 0*174000 03 SQ.FT WEIGHT * 0.265000 04 LBS 


STATIC PERFORMANCE AT AN ALTITUDE - 0.0 FT 

WITH MINIMUM TIME AND MOST ECONOMICAL CLIMB SCHEDULES TO A FINAL ALTITUDE • 0.100000 05 FT 


MINIMUM LEVEL FLIGHT SPEED - 0.904650 02 FTSSEC 
LIFT COEFFICIENT • 0.156380 01 ORAG COEFFICIENT • 0.484180 00 

MAXIMUM LEVEL FLIGHT SPEED » 0.252570 03 FT/SEC 
LIFT COEFFICIENT - 0.200620 00 ORAG COEFFICIENT • 0.290630-01 


MAXIMUM CLIMB ANGLE - 0.10220D 02 0E6 
VELOCITY FOR MAXIMUM CLIMB ANGLE « 0.117100 03 FT/SEC 
LIFT COEFFICIENT • C. 933380 00 DRAG COEFFICIENT « 0.854750-01 


VELOCITY FOR MAXIMUM ENDURANCE - 0.125710 03 FT/SEC 
POWER FOR MAXIMUM ENDURANCE • 0.275450 05 FT-LBS/SEC 
LIFT COEFFICIENT * 0.809850 00 DRAG COEFFICIENT - 0.669610-01 


VELOCITY FOR CLASSICAL MAXIMUM RANGE - 0.142050 03 FT/SEC 
LIFT COEFFICIENT - 0.634230 00 DRAG COEFFICIENT - 0.496190-01 


SERVICE CEILING - 0.194420 05 FT 
VELOCITY AT SERVICE CEILING - 0.175540 03 FT/SEC 
LIFT COEFFICIENT - 0.764240 00 ORAG COEFFICIENT • 0.616530-01 

ABSOLUTE CEILING - 0.212360 05 FT 
VELOCITY AT ABSOLUTE CEILING * 0.18023D 03 FT/SEC 
LIFT COEFFICIENT - 0.770560 00 DRAG COEFFICIENT - 0.623480-01 


MAXIMUM RATE OF CLIMB SCHEDULE FROM 0.0 FT TC 0.10000D 05 FT 


HIFTI 


R/CIFT/SECI 

VfFT/SECt 

PIFT-LBS/SEC) 

CL 


CO 

TISEC1 


0.0 


0.222570 

02 

0.136010 

03 

0.673540 

05 

0.691830 

00 

0.544600-01 

0.0 


0.500000 

03 

0.216590 

02 

0.136620 

03 

0.859210 

05 

0.695810 

00 

0.548210-01 

0.2277SD 

02 

0.100000 

04 

0.210650 

02 

0.13725D 

03 

0.845040 

05 

0.699660 

00 

0.551740-01 

0.661850 

02 

0.150000 

04 

0.204750 

02 

0.137890 

03 

0.83102D 

05 

0.703390 

00 

0.555190-01 

0.702630 

02 

0.200000 

04 

0.198890 

02 

0.13856D 

03 

0.117160 

05 

0.707000 

00 

0.558560-01 

0.950 4 20 

02 

0.250000 

04 

0.193070 

02 

0.139250 

03 

0.803450 

05 

0.710490 

00 

0.561S5D-01 

0.120560 

03 

0.300000 

04 

0.187280 

02 

0.139960 

03 

0.789890 

05 

0.713840 

00 

0.565040-01 

0.146860 

03 

0.350000 

04 

0.181540 

02 

0.140690 

03 

0.77 6 4 70 

05 

0.717080 

00 

0.568 140-01 

0.173980 

03 

0.40000D 

04 

0.175830 

02 

0.141450 

03 

0.763200 

05 

0.720180 

00 

0.571150-01 

0.201970 

03 

0.450000 

04 

0.170150 

02 

0.142220 

03 

0.73C080 

05 

0.723160 

00 

0.57405D-01 

0. 230 8 80 

03 

0.500000 

04 

0.164510 

02 

0.143020 

03 

0.737090 

05 

0.726000 

00 

0.576840-01 

0.260770 

03 

0.550000 

04 

0.158910 

02 

0.143840 

03 

0.724250 

05 

0.728720 

00 

0.579530-01 

0.291700 

03 

0.600000 

04 

0.153350 

02 

0.144680 

03 

0.711550 

05 

0.731310 

00 

0.582110-01 

0.323730 

03 

0.650000 

04 

0.147820 

02 

0.145540 

03 

0.698990 

05 

Q. 733760 

00 

0.584570-01 

0.336950 

03 

0.700000 

04 

0.142330 

02 

0.146430 

03 

0.686570 

05 

0.736080 

00 

0.586910-01 

0.391420 

03 

0.750000 

04 

0.136870 

02 

0.147340 

03 

0.674280 

05 

0.738280 

00 

0. 589140-01 

0.427260 

03 

0.800000 

04 

0.131440 

02 

0.148280 

03 

0.662130 

05 

0.T40340 

00 

0.591240-01 

0.46454D 

03 

0.B5000D 

04 

0.126060 

02 

0.149240 

03 

0.650110 

05 

0.742270 

00 

0.593220-01 

0.503390 

03 

0.900000 

04 

0.120700 

02 

0.150220 

03 

0.638230 

05 

0.74407D 

00 

0.59508D-01 

0.543940 

03 

0.950000 

04 

0.115390 

02 

0.151220 

03 

0.626490 

05 

0.745730 

00 

0.596810-01 

0.586320 

03 

0.100000 

05 

0.110100 

02 

0.152260 

03 

0.614880 

05 

0.747270 

00 

0.598410-01 

0.630690 

03 


Table 2. Performance of the Cessna 182 with the 
the general drag polar. 
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P ( sea level ) 



Figure 3. Power required and power available versus velocity. 






2000 2200 2400 2600 2800 3000* 

Weight (pounds) 



Altitude (feet) 


Figure 4. Effects of weight and altitude on 
maximum and minimum level flight 
speeds and maximum rate of climb. 


Maximum Rate of Climb (ft. /sec) 


Weight ( pounds ) 


Figure 5 


Service and absolute ceilings 
for various weights. 




POWER AVAILABLE VS. VELOCITY 
REFERENCE ALTITUDE - 0.0 FEET 

PAI FT-LB S/SEC I VIFT/SECI 

O.C 0.0 

0.363000 05 0.281500 02 

0.653400 05 0.563000 02 

0.671200 05 0.84450D 02 

0. 101640 06 0.112600 03 

0 .106500 06 0.140750 03 

0.113260 06 0.168900 03 

0.117610 06 0.197050 03 

0.119790 06 0.225200 03 

0.120520 06 0.253350 03 

0. 12 1240 06 0.281500 03 

0.121970 06 0.309650 03 

0.123420 06 0.337800 03 

0.123420 06 0.36595D 03 

0.123420 06 0.394100 03 


AIRCRAFT CHARACTERISTICS 

CD - 0.470000-01 ♦ 0. 10153D-01*CL**2 ♦ 0.434780-014CL** 0.541160 01 
WING AREA - 0.160000 03 SQ.FT WEIGHT - 0.275000 04 LBS 


STATIC PERFORMANCE AT AN ALTITUDE - 0.0 FT 

WITH MINIMUM TIME ANO MOST ECONOMICAL CLIMB SCHEDULES TO A FINAL ALTITUDE - 0.100000 05 FT 

MINIMUM LEVEL FLIGHT SPEED * 0.90382D 02 FT/SEC 
LIFT COEFFICIENT - C. 157160 01 DRAG COEFFICIENT - 0.574190 00 

MAXIMUM LEVEL FLIGHT SPEED - 0.227310 03 FT/SEC 
LIFT COEFFICIENT - C. 246460 00 DRAG COEFFICIENT « 0.47650D-01 


MAXIMUM CLIMB ANGLE - 0.130610 02 DEG 
VELOCITY FOR MAXIMUM CLIMB ANGLE - 0.117270 03 FT/SEC 
LIFT COEFFICIENT - 0.933490 00 DRAG COEFFICIENT ■ 0.858060-01 


VELOCITY FOR MAXIMUM ENDURANCE - 0.123510 03 FT/SEC 
POWER FOR MAXIMUM ENDURANCE - 0.287720 05 F7-LBS/SEC 
t I F7 COEFFICIENT • C.B41640 00 DRAG COEFFICIENT - 0.712960-01 


VELOCITY FOR CLASSICAL MAXIMUM RANGE » 0.130600 03 FT/SEC 
LIFT COEFFICIENT - C. 752760 00 DRAG COEFFICIENT - 0.62I03C-01 


SERVICE CEILING - 0.221060 05 FT 
VELOCITY AT SERVICE CEILING - 0.177820 03 FT/SEC 
LIFT COEFFICIENT - O.81610D 00 DRAG COEFFICIENT * 0.68464C-01 

ABSOLUTE CEILING ■ 0.237250 05 FT 
VELOCITY AT ABSOLUTE CEILING - 0.182620 03 FT/SEC 
LIFT COEFFICIENT - C.82C440 00 ORAG COEFFICIENT - 0.687320-01 


MAXIMUM RATE OF CLIMB SCHEOULE FROM 0.0 FT TO 0.10000D 05 FT 


HIFTI 


R/CIFT/SEC) 

V( FT/SEC 1 

PCFT-LBS/SEC) 

CL 


CO 

T(SEC) 


0.0 


0.28095D 

02 

0.130500 

03 

0.106870 

06 

0.75389D 

00 

0.621960-01 

0.0 


0.500000 

03 

0.274160 

02 

0.131160 

03 

0.105150 

06 

0.757270 

00 

0 .624790-01 

0.180170 

02 

0. 100000 

04 

0.267410 

02 

0.131840 

03 

0.103460 

06 

0.760550 

00 

0. 627580-01 

0.364850 

02 

0.15000D 

04 

0.26070D 

02 

0.13254D 

03 

0.101770 

06 

0.763710 

00 

0.630310-01 

0.554240 

02 

0.20000D 

04 

0.254020 

02 

0.133260 

03 

0.100110 

06 

0.766760 

00 

0. 632990-01 

0.748550 

02 

0.25000D 

04 

0.24739C 

02 

0.134000 

03 

0.984620 

05 

0.769700 

00 

0.635610-01 

0 .948020 

02 

0.300000 

04 

0.240790 

02 

0.134750 

03 

0.968290 

05 

0.772530 

00 

0.638170-01 

0.115290 

03 

0.350000 

04 

0.234230 

02 

0.135320 

03 

0.952130 

05 

0.775250 

00 

0.640660-01 

0.136350 

03 

0.400000 

04 

0.227710 

02 

0.136320 

03 

0.936120 

05 

0.777860 

00 

0.64308D— 01 

0.158000 

03 

0.450000 

04 

0.221230 

02 

0.137120 

03 

0.920270 

05 

0.780360 

00 

0.645440-01 

0.180280 

03 

0.500000 

04 

0.214780 

02 

0.137950 

03 

0.904570 

05 

0. 782750 

00 

0.647710-01 

0.203220 

03 

0.550000 

04 

0.208370 

02 

0.138800 

03 

0.889030 

05 

0.78503D 

00 

0.649910-01 

0.226850 

03 

0.600000 

04 

0.202000 

02 

0.139670 

03 

0.873640 

05 

0.7B721D 

00 

0.65203D-01 

0.251230 

03 

0.650000 

04 

0.195670 

02 

0.140550 

03 

0.858400 

05 

0.789280 

00 

0.65407D— 01 

0.276380 

03 

0.700000 

04 

0.189380 

02 

0.141460 

03 

0.843310 

05 

0.791250 

00 

0. 656020-01 

0.302360 

03 

0.750000 

04 

0.183120 

02 

0.142380 

03 

0.828360 

05 

0.793130 

00 

0.657910-01 

0.329210 

03 

0.000000 

04 

0.176900 

02 

0.143320 

03 

0.813560 

05 

0.79493D 

00 

0. 659730-01 

0.357000 

03 

0.850000 

04 

0.170710 

02 

0.144280 

03 

0.798910 

05 

0.796650 

00 

0.66 149D-0 1 

0.385770 

03 

0.900000 

04 

0.16457D 

02 

0.143250 

03 

0.784410 

05 

0.798280 

00 

0.66317D— 01 

0.415610 

03 

0.950000 

04 

0.158460 

02 

0.146250 

03 

0.770040 

05 

0.799830 

00 

0.664780-01 

0.446580 

03 

0.100000 

05 

0.152390 

02 

0.147260 

03 

0.755830 

05 

0.801310 

00 

0.666310-01 

0.476760 

03 


Table 3. Performance of the Navion for a general drag polar. 
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TAKE-OFF AND LANDING PERFORMANCE 


The static performance previously discussed has been concerned primarily 
with the cruising and climbing aircraft. However, the complete performance 
analysis must also include a discussion of take-off and landing. Take-off 
performance analysis can usually be divided into two parts, ground run and 
climb over a 50 foot obstacle. Analogously, landing consists of the approach 
(from a 50 foot altitude to touch-down) and the ground run. A brief analysis 
of both take-off and landing performance taken largely from Reference 15 is 
presented below. 



Figure 6. Forces acting on an aircraft during landing and take-off. 


Take-Off Ground Run 


The distance traversed by an airplane on level ground in accelerating 
from one speed to another can be expressed by 


where 



> 



and a = 


dV_ 

dt 


The ground run can thus be found by integrating the above expression to 
yield 


S 


G ~ 




( 21 ) 


If the airplane is starting from rest V x = 0 and V y = V^qf = velocity at 
lift-off. The effects of wind on take-off can be accounted for by letting 
V x = V w = velocity of the headwind and by replacing V by (V - V w ) since a 
is a function of ground velocity (V - V w ). 
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The frictional force of the wheels on the runway is a force which acts 
in the same direction as the drag force. This force is proportional to the 
weight less the lift and is given by 

F f = (W - L) vi 

where p is the coefficient of friction. The value of p depends upon the type 
of runway surface used. Typical values of p are-given below (Ref. 29): 


SURFACE 

y 

Concrete, asphalt, or wood 

0.02 

Hard turf 

0.04 

Average field - short grass 

0.05 

Average field - long grass 

0.10 

Soft ground 

0. 10 -»■ 0.30 


Table 4. Typical values of p for various runway surfaces. 


In addition to the frictional force opposing the thrust there is a drag 
force so that the acceleration during the ground run is given by 


_ g[T - D - P(W - L)3 
a _ _ 


( 22 ) 


Expressing drag 
be expressed as 


and lift in coefficient form the distance equation 



w (V - V w ) d V 

? [T - pW - (C D - pC L ) l p 0 aV 2 ] 


9 


(21) can 


(23) 


where Sq is the ground distance required for take-off. For take-off Cq, C|_, 
and W can be assumed to be independent of velocity and the thrust, T, can 
be found from P/V where P wi I I be the maximum power. Since dt = dV/a", the 
time to li*+-off can be expressed as: 



W dV 

g [T - pW - (C D - pC L ) i p 0 aV 2 s3 


(24) 
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Since a numerical integration technique is required to evalaute the 
integrals given in Equations (23) and (24), a comparatively simple method 
for predicting the take-off ground run and the time to lift-off is also 
offered. This method (Ref. 14), presented in 1933, was intended as a 
practical approximation to a difficult problem. The steps proceed as 
fo I lows: 


(1) Evaluate 


( k T )(BHP) 
o 


and 


W (W) (N) (D) 

JF _ (550)(thpm) ( th P _ £ 
W W V s thp m " L 


where: 


N 

D 

W 

V 

+hp m 


static thrust coefficient (Figure 7 or 8) 

engine brake horsepower 

engine speed in revolutions per minute 

propeller diameter in feet 

take-off weight in pounds 

coefficient of wheel friction (Table 4) 

maximum thrust horsepower 


) = ratio of thrust horsepower at speed V 
™Pm to thrust horsepower at maximum speed 
(Figure 9) 


= the reciprocal of the maximum value of 
L/D which is either known or estimated 
(Figure 10) 

V s = take-off velocity. 


(2) Evaluate 


c - K s v a 

° T i 
(— ) 
w 


and 


*o ~ 


Vs 

(li, 

w 


where: 

S Q = take-off ground run 
t 0 = time required to lift-off 
K s = take-off coefficient in Figure 1 1 as a 
function of (Tp/W)/(T^/W) 

Kj. = time coefficient in Figure 12 as a function 
of (T F /W)/nVW). 


(3) 


Corrections for take-off headwind and take-off weight changes 
can be obtained from Figures 13 and 14 by using the two relations 
given below: 


S w = S o 


and 



W1 

w 


2 
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where: 


s w 

5 o 


F = 


ratio of take-off distance in headwind to 
distance in a calm (Figure 13) 
weight correction factor (Figure 14). 


Because of the lack of good theoretical lift, drag, and thrust data a 
comparison of the two methods presented for estimating the take-off ground 
distance and the time to lift-off has been omitted. 



Figure 7. Static thrust coeffi- 
cient versus 0. 



Figure 8. Static thrust coeffi- 
cient versus V/nD. 
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F ino! ne t thrus t 
,0 ‘ Initial net thrust 


* Maximum speed V„ 

t 

.02 

— i — i — i i — i i i i i i 

.03 .04 .05 .06 .07 

Parasite droq coefficient t C Dp 

Figure 9. General full-throttle 



t.hp. curve. 

Figure 

10. Maximum L/D versus 


parasite drag 
coeff icient. 


111 


■ill 

mmm 


0/6 .030 .034 .038 .032 036 .040 

Tohe~off coe fftcient, K, 


Figure 11. Take-off run in calm. 


Final net thrust w L 


* Initial net thrust l T 0 


Figure 12. Coefficient for time 
to take-off. 
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Climb to 50 Feet 


The approximate ground distance traveled in attaining an altitude of 
50 feet after lift-off can be expressed by a relation similar to the one 
given for the ground run distance: 




a 


dt 


P 


or 


>50 



w (v - V dV 

g [T - (C D ) i p 0 av 2 ] 


(25) 


where V 5 q is the velocity of the airplane at an altitude of 50 feet. The 
time required to attain an altitude of 50 feet can similarly be expressed 
as: 


t = f V5 ° W dV 

50 ' J g ET - CC D > i p av's] • 

V/ u o 


(26) 


The total take-off distance traversed in going from a position of rest to an 
altitude of 50 feet is thus Sq + S^q, and the total take-off time is tg + t^Q 
Equations (3), (4), (5), and ( 6 ) can be integrated numerically if the veloc- 
ities, which are the limits of integration, are known; if Cq and Cl are 
considered independent of velocity, and if the functional form of T(V) is 
known. As an aid for estimating V^qp and V 5 Q it should be noted that Part 23 
of the Federal Aviation Regulations (Ref. 30) requires that V 50 t> e at least 
1.3 times the zero thrust stall speed. One would also expect V^Qp to be 
approximately 1.1 times the stall speed. 


Landing Approach 

The landing approach can be divided into two parts— -a steady-state glide 
path where the airplane is in the final landing configuration prior to 
touch-down and the flare. 
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Touchdown 




Figure 15. Typical pattern for landing approach, flare, and touchdown. 


Assuming that the flare is a circular arc (Figure 15) and that y is small, 
then S^, the ground distance traversed in descending from an altitude of 50 
feet to touchdown, can be expressed as: 



(27) 


Since y is small, V is assumed to be constant during the approach glide. An 
expression for y and R must now be determined. With the weight approximately 
equal to the lift, y can be written as 


y 


Cd 

c L 


I 

w 


(28) 


The acceleration normal to the flight path needed to flare is attained by 
rotating the airplane to a higher C|_ value, say C|_. The lift force during 
the flare can be written 
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L' = c’p o | SV 2 . 

Then the force normal to the flight path is 

= l' - W cos y = L* - W . 

However, during the glide L - W; thus l'/W = c'/C . The normal force can 
therefore be written as L L 

c * 

F N = w ^ - W = W(n - 1) (29) 

where n = C./C. = load factor (the maximum value of n which may be applied to 

C[_ is dictated by stall or buffet limits). 

The force normal to the flight path can also be expressed as 


F N = £ a N 


where 


Equating the two expression for normal force yields the relation 

W(n - 1) = - ~ 

g R 


g(n - 1) 

BU+ ’ V 2 - nW _ 2 

S Po aC- S p Q aC’ ‘ 

With equation (28), (30), and (31) and the fact that C^_ = nC^, equation (27) 
can be expressed as 

W , C D T. 

, 80 . ^ ( CT " r 

S = ~2il + L . (32 

( C D Ts p Q ag(n - DC. 

c L - w 

The ground distance for approach can thus be estimated if Cq, C l , T, and 
the change in lift coefficient which occurs during the flare, n, are known. 


Landing Ground Run 


The landing ground run can be described as (1) a short ground run 
(approximately two seconds) immediately following touchdown while the airplane 
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is being changed from landing configuration to braking configuration and 
(2) the remaining ground run which brings the airplane to a complete stop. 
The distance covered in landing transition is taken to be 


S TRAN " 


V TD + V B 

2 At TRAN 


(33) 


where 

Vjd = touchdown velocity* ** 

Vg = speed at full braking configuration, 

At TRA M = +he transition time from touchdown to full 
braking (approximately 2 seconds). 


The ground run braking distance is obtained in the same manner as the take- 
off ground run, 


S 


B = 



d V 


(34) 


where 



a 


= initial braking speed, 
= headwind velocity, 

= deceleration rate. 


The acceleration (deceleration if it is negative) term derived in the take-off 
section may be used here with an appropriate braking value of y (a good 
approximate value of y for aircraft tires on an asphalt runway is 0.3) 

Thus, 


• V w W(V - V )dV 

w 

v gCT - y B w - (C D - y B c L ) ^ V 2 p 0 s] 
v B 


(35) 


In using the above equation the values of Cq and C|_ must be those actually 
obtained in a landing run. For instance, if the flaps are retracted after 
touchdown then the approach C L should not be used to calculate the braking 
distance. Similarly if an aircraft is equipped with reversible-pitch 
propel lers, the thrust term should take this into account. If a reversible 
propeller is not used then T should represent the idle thrust. 

If one desires to achieve the minimum braking distance then a large 
negative value of the acceleration term is desired. This is best obtained 
with a negative thrust (reversible propeller), a high drag coefficient, and 
a low value of Cj_. 


* It should be remembered that the touchdown velocity must be low enough 

that the lift will not be greater than the weight at the touchdown altitude. 
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PATH PERFORMANCE 


DESCRIPTION OF THE PATH PERFORMANCE EQUATIONS AND THE INTEGRATION PROCEDURE 
EMPLOYED TO OBTAIN FLIGHT TIME HISTORIES 


In recent years, development of the digital computer and improved numeri- 
cal techniques have made possible the prediction of aircraft performance in 
a somewhat more sophisticated manner. Prior to this, the static analysis 
detailed in the preceding section was used exclusively in estimating aircraft 
performance. However, when more-refined analysis is desired the approach 
most often employed is integration of the vehicle equations of motion, a 
system of simultaneous, ordinary differential equations. This integration 
yields time histories of both vehicle and flight path parameters for various 
control inputs. By repeatedly solving these equations for diversified inputs, 
the practical optimum of some desired parameter ( i.e ., range, rate of climb, . .) 
may be approached. 

The general equations governing aircraft performance as derived in 
Appendix B have been programmed for numerical solution. A discussion of the 
solution procedure and of some typical results for various input combinations 
constitutes the majority of the section. A detailed user guide, complete 
with an example case and a program listing, is presented in Appendix D. 

The basic equations (B-19) are presented below for reference. 


x = V cos y (36) 

• 

h = V sin y (37) 

• gP g SV 2 

V = WV - W — C D (a) p(h) - g sin y (38) 

y = J y- c L (a) p(h) “ y - cos y (39) 

W = - cP (40) 

p(h) = p 0 ( 1 . 0 - 6.86 x 10“ 6 h) 4 * 26 (41) 

C, (a) = C, a + C L (42) 

a (a=0) 

, k 7 

C D (a) = C Dq + k[C L (a)] 2 + k 1 [C L (a)] 2 (43) 

P < P max (h,V) (44) 
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The first equation, which relates the horizontal distance traveled to the 
velocity and flight path angle, is added to those of Appendix B to permit 
the direct calculation of range. Note that the fuel-flow rate is considered 
to be directly proportional to engine power rather than to thrust, since 
concern here is primarily for aircraft ppwered by piston engines. Also, 
the variation in thrust with angle of attack is restricted to small a, since 
the thrust vector is assumed to lie along the body axis.* 

In programming the general equations, Equation (43) was simplified by 
setting k equal to zero, yielding the following functional form. 

c D ( a) = C Dq + k 1 [C 1 _(a)] k 2 (45) 

When drag data are available, the three parameters in Equation (45) (Cq , 
k^, k 2 > may be found by the curve-fitting scheme presented in the previSus 
section. A comparison of the accuracy obtained by fitting drag data with 
the three parameters of Equation (45) as opposed to the four parameters of 
Equation (43) indicates a negligible difference except for angles of attack 
approaching stall. Since most performance trajectories do not entail con- 
tinued operation near the minimum speed, this simplification introduces no 
significant error. For calculations in which operation near the minimum 
velocity is of primary interest, the reader should investigate the procedure 
presented in the previous section which employs the four parameter fit to 
drag data. 

If drag data are not accessible, the standard parabolic drag polar 
evolves from Equation (45) as fol lows: 


where 


C D (a) " C D 0 



+ — — 
TreAR 


e = Oswald's efficiency factor 
AR = wing aspect ratio 
Cq = drag coefficient which is independent 
of angle of attack. 


(46) 


From a mathematical viewpoint, the performance problem is one of solving 
four nonlinear ordinary differential equations involving six dependent 
variables and one independent variable — time. Consequently, there are two 
degrees of freedom which imply an infinite number of trajectories for each 
set of initial conditions. To obtain a unique trajectory, time histories 
of six unknowns must be specified. Although there are fifteen possible 
combinations of these unknowns taken two at a time, the programmed solution 
procedure allows the user to specify only fourteen sets of inputs since 
weight and power can not be specified independently. The optimal desired 


* These restrictions are not fundamental to the technique. More general 
expressions, if available, may easily be substituted without significant 
increase in computing time. 
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(range,- rate of climb, etc.) strongly influences the selection of variables 
to be specified. The variation of any chosen variable and its derivatives 
is assumed known throughout the entire trajectory. Initial conditions for 
each of the differentiated variables must be provided for use in starting 
the integration procedure. 

The integration technique, essentially a Runge-Kutta Adams-Bashforth 
predictor-corrector procedure, was altered on the basis of a paper by Charles 
Treanor (Ref. 31). Since it is self-starting, the modified Runge-Kutta was 
used to calculate the first four points. Then the modified predictor- 
corrector performed the remaining integration until the step size was either 
increased or decreased due to some error criterion. 


The following is a brief synopsis of the evolution of the integration 
procedure from the classical Runge-Kutta into its final form. The modifica- 
tion proposed by Treanor is applicable to any differential equation in which 
the derivative of the dependent variable has a strong dependence on the 
difference between the value of the variable and that of some slowly changing 
function. This type behavior is notably evident in calculating the derivative 
of flight path angle as can be seen by grouping the right hand side of 
Equation (39) as fol lows: 


; = £ ( £<h)VfS } !_ £ 

T V 2 L W v Y 


(47) 


Substituting for the lift and rearranging yields: 

Y = - f [cos Y " . (48) 

Obviously, the derivative is strongly dependent upon y and the slowly 
changing lift to weight ratio. When solving a system of equations using 
standard Runge-Kutta Adams-Bashforth technique, the above equation induces 
enormous oscillations which severely restrict the size of the integration 
increment. The Treanor modification to the classical Runge-Kutta method 
removes these oscillations in any integration interval where this strong 
dependence is detected. For regions not exhibiting this particular 
behavior, the modified method reverts to classical fourth order Runge-Kutta. 

However, even with the enlarged integration step size permitted by the 
modified Runge-Kutta, solution of the system of equations proved to be very 
costly (in computer time) when integrating trajectories of long duration 
(i.e., eight hours or more). This is because Runge-Kutta procedures 
essentially obtain the solution twice for each integration interval. In 
a particular interval the results are first computed at t 2 = + At and 

then recomputed using steps of t/2. A comparison between the two results 
indicates the accuracy obtained and determines whether the integration step 
size should be increased or decreased. Thus the equations are actually 
solved twice, once in normal step sizes and again in half steps with a 
comparison following the larger step. This computational complexity 
suggests the use of a fourth order predictor-corrector such as that of 
Adams-Bashforth (Ref. 32). A modification of this predictor-corrector. 
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which can accommodate derivatives strongly related to the value of the depen- 
dent variable and which provides the compatibility conditions necessary for 
use with the modified Runge-Kutta procedure, is derived in Appendix G. Use 
of this combined integration procedure permitted the step size to increase 
by one to three orders of magnitude depending upon the type and nature of 
the specified variables. Computational times were so drastically reduced 
that integration of eight to ten hour flight trajectories were relatively 
inexpensive. Execution times for an eight hour trajectory may vary from 
about nine seconds, for a solution having altitude and angle of attack 
specified as constant, to slightly more than two minutes for a trajectory 
with a superimposed oscillatory mode which results from specification of 
angle of attack as a constant and power as the maximum available. 

The constraint on the general performance problem. Equation (44), demands 
that the power for any flight maneuver always be less than or equal to the 
maximum power avai lable at a particular altitude and velocity. To fulfi I I 
this requirement, the maximum power available as a function of altitude and 
velocity must be determined for the aircraft under investigation. The 
calculation of P max (h>V) for use in this general integration procedure is 
detailed in the previous section. Briefly, the technique utilizes a spline 
fit of several power available versus velocity points at some reference 
altitude to predict the maximum power available at any altitude and velocity. 
The program compares the actual power calculated with P max (h,V) at each 
integration step to insure that the resulting trajectory satisfies the 
maximum power available constraint. 


I 
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TYPICAL RESULTS FROM INTEGRATION OF THE PATH PERFORMANCE EQUATIONS 


The performance of a vehicle in actual flight may exhibit many features 
of interest in addition to the equilibrium maximums predicted by static 
calculations. Some aircraft experience accelerations too severe for reliable 
static approximation. For example, in a zoom maneuver or a rapid climb, 
the airframe dynamics are obviously nonlinear and their description demands 
a more sophisticated treatment. Large weight changes associated with high 
power loadings also emphasize the need for this type analysis. The following 
examples serve to illustrate this treatment and will, perhaps, encourage the 
reader's use of the more realistic path performance approach possibly 
employing the accompanying computer program (Appendix D). 

The following examples all employ the classical parabolic drag polar. 
However, computation with the three parameter polar of Equation (45) proceeds 
with equal ease. A value for Cd q was taken from Reference 1 and c, the 
specific fuel comsumption, was given an average value of 0.6 Ibs/hp-hr. 

First, the determination of aircraft range for several different sets 
of specified parameters was investigated. When considering the prediction 
of maximum range for a given amount of fuel, the investigator who draws upon 
his previous experience immediately considers angle of attack as one of the 
specified variables, since for quasi-steady conditions flying at a = a (L/D) max 
produces the greatest range per pound of fuel. Calculations carried out 
using the path performance equations (and discussed below) prove the validity 
of this result for light aircraft. The second specified variable could be 
any of several whose effect on the resulting trajectory might yield an in- 
creased range. To determine a(L/D) max for The +h ree parameter drag polar 
of Equation (45), first form the raiio (C|_/C D ) and differentiate with respect 
to C L , set the derivative equal to zero, then solve for C|_ ^ l/q ) • This 

process yields max 

C L(L/D) max 

For the parabolic drag polar, \<2 = 2, and 

0 -m’ /2 

l - ( L/I»ma X U, / 




Inserting (50) into Equation (42) yields a for maximum lift to drag ratio. 
Each trajectory began with similar initial conditions and was terminated 
upon consumption of 222 pounds of fuel. Initial and final values for 
several pairs of specified variables are given in Table 5. 
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Table 5. Results of several test cases designed to obtain maximum range. 



The specified variables for Case (1) are constant angle of attack and 
constant altitude. During the eight-hour trajectory, the weight decreased 
by 222 pounds with an associated velocity drop of 6.4 feet per second and 
a corresponding power reduction of about 6 hp. In Case (2) a smaller angle 
of attack (a = 3 degrees) was specified and the other parameters remained 
as in Case (1). After expending 222 pounds of fuel, the range was 43 miles 
less than for Case (1). This illustrates the advantage of flying with angle 
of attack for maximum lift to drag ratio. In Case (3), angle of attack and 
a siightiy larger velocity than that resulting from previous trajectories 
were specified. This attempt to further range through a slight increase in 
speed proved futile, since the vehicle climbed nearly 2,000 feet and then 
began to glide after using the allotted fuel until at the same final altitude 
it obtained an equivalent range. 

In Case (4) the power was specified as the minimum required for level 
flight — at 10,000 feet with a = ot( L/D) ma x — mu I i - i p I i ed by a factor to allow 
for density variation as the vehicle climbs due to fuel usage. Once the 
weight reached 2,428 pounds (222 pounds of fuel burned), the power was set 
equal zero and the vehicle permitted to glide until it descended to 10,000 
feet altitude. A similar power schedule was followed in Case (5) with the 
exception that power was assigned to be the maximum available as a function 
of local altitude and velocity. As before, when fuel consumption reached 
222 pounds, the power was set equal zero with the effect being a glide 
descent to 10,000 feet. Figure 16 presents a comparison of power schedules 
with the resulting altitude and velocity variations. The vehicle of Case (4) 
glided the last six miles; whereas, the Case (5) craft finished by gliding 
for nearly 45 miles. However, as Table 5 indicates, both trajectories 
produced the same range, but Case (5) required nearly 1.5 hours less flight 
time. 


Motivation produced by the performance optimization study of Reference 33 
led to the investigation of aircraft range using a combination of three 
power settings. The trajectory was divided into three sections — ascent, 
cruise, and descent — with a different power setting corresponding to each 
regime. First, the power was designated as the maximum available for the 
ascent stage; then at some predetermined altitude the power was reset so 
that during cruise the thrust and drag were equal; finally, when the assigned 
portion of fuel (common to each trajectory) was consumed, the power was shut 
off and the vehicle executed a glide descent back to its initial altitude. 
Since range was of primary concern, the second variable was specified as 
a = a (L/D)max' A family of trajectories was integrated with the altitude 
for switching the power from maximum to that for thrust equal drag being 
varied from 14,000 feet to 22,000 feet in increments of 2000 feet. Figure 
17 presents altitude and power variations corresponding to the upper and 
lower switching altitudes. After descent to 10,000 feet, the range for 
every trajectory was equivalent. Since these flights were only five hours 
in length, the final range was compared with Case (1) of Table 5 at the 
five hour mark and they were also equivalent. The major conclusion from 
the above investigation of aircraft range is that most trajectories will 
produce ranges of equal magnitude whenever one of the specified variables 
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Horsepower 
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Figure 16. Time histories of several parameters resulting from analysis of maximum range. 


NOTE 5 Power switched from P mox to that for thrust 
equal drag at an altitude of : 

14000 ft. 



Time histories of power and altitude for trajectories using 
three different power settings. 



is set as a = cx.( l/D ) max* Thus the major objective for most pilots desiring 
maximum range should be flight at approximately the angle of attack for 
maximum lift to drag ratio. To obtain more realistic results than those in 
Table 5 (i.e. 795 miles on 222 pounds of fuel is quite an exaggeration) 
and to accurately approximate a(|_/D) max > improvement in the prediction of 
three-dimensional lift and drag characteristics is mandatory. More precise 
knowledge of the variation of fuel flow rate with power required, aircraft 
speed, and altitude, and engine speed and throttle setting is also necessary. 
Acceptable trends and desired objectives may be formulated using the classical 
lift and drag relations; however, realistic values from performance calcula- 
tions will result only when high quality aerodynamic and engine test data 
or sophisticated prediction techniques are utilized to estimate the body 
forces and the fuel flow rate. 

Since the objective of the preceding analysis was determination of 
parameter behavior throughout the flight, large time scales of several hours 
or more were employed. These large time divisions, however, often mask 
some interesting behavior which occurs simultaneously but can be observed 
only on a small time scale. An inset, describing the first two minutes of 
flight path angle variation with time for the trajectory of Case (5), is 
presented in Figure 18. 
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Figure 18. Sample phugoid trajectory. 




* This result should not be too surprising. Consider for the moment the 
situation which would exist if the drag were zero. Then the energy required 
to traverse any trajectory between point A and point B would be the same, 
namely zero. Because'of the existence of drag, energy can be dissipated in 
two ways: that required to overcome friction in steady flight and that 

associated with the increasing disorder (entropy) of accelerated motion. Now, 
if the drag is relatively small to begin with and the accelerations are also 
modest, then it is reasonable to expect that a great many different flight 
trajectories will have very similar energy requirements, that is, they will 
require the same amount of fuel plus or minus a pound or two. One would expect 
very different results for trajectories involving near sonic speeds with 
significant accelerations. 
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Since this behavior normally results whenever angle of attack (i.e. lift 
coefficient) is held constant, the oscillation observed in Figure 18 is 
commonly known as a phugoid trajectory. The period of oscillation, about 
21 seconds, is typicaf of the phugoid motion for light aircraft. The 
oscillations are damped as the velocity increases and the vehicle climbs, 
due to the decrease in amount of excess power available. Observation of 
this trajectory on a larger time scale, reveals only an average y which 
decreases from about 4.4 degrees to approximately zero late in the flight. 
Recovery of the aforementioned phugoid motions from the total trajectory 
symbolizes the very general nature of this solution and emphasizes the 
lack of most restrictive assumptions. 

The application of a general solution technique to landing flight 
appears very informative. For analysis, the landing in subdivided into two 
portions, the approach and the flare, with velocity and flight path angle 
specified in each interval. A basic geometrical description of the glide 
path and flare is presented in Figure 19. 



Figure 19. Geometrical description of landing analysis. 


The overall concept is first to descend at constant velocity and flight path 
angle until reaching an altitude of and then execute a flare with velocity 
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decreasing and y becoming zero at touchdown. For illustration, the flight 
path angle is specified so as to produce a circular arc flare and velocity 
is linearly decreased from the approach speed to some touchdown velocity. 
Specification of y in the flare region necessitates finding the equation of 
a circle of radius R with center at (j, k) passing through points (&j, h^) 
and (Hj t &2 + b » 0) • This equation is of the form: 


Since 


(x - j) 2 + (h - k) 2 = R 2 . 

(52) 

tan Y = 

dx (h - k) 

(53) 

flare regime is given by: 


Y = - tan 1 (jpE-^) . 

(54) 


For calculation of j and k, values of h Q (altitude at beginning of approach), 
y 0 (constant glide path angle during approach), and either h^ (altitude of 
flare initiation) or b (distance from intersection of extended approach path 
with ground to the point of touchdown) must be given. Using these three 
parameters, values for j and k may be determined from the following formulae: 


h| = b sin y 0 (55) 

R = b/tan(y c /2) (56) 

Jbj = ( ^o “ hj)/tan Yo (57) 

k = R (58) 

j = JL + /hj (2R - h t ) . (59) 


Thus velocity and flight path angle may be specified as follows: 


Y = 


" Y, 



hi < h < h Q 


0 £ h £ h 1 


(60) 


V = 


V 


o 


<V o - V T>‘^) + V T 


h 1 < h « h o 


0 ^ h < hj 


(61) 
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where 


Also, 


V Q = approach velocity 
Vj = touchdown velocity. 


0.0 




h. < h < h 
1 o 

0 < h < h 1 


( 62 ) 


h. < h < h 
l o 

0 £ h < h 1 

Consider the following example for which: 
h 0 = 1500 feet 

Y q = 2.5 degrees = 0.04363 radians 
b = 400 feet 
V Q = 140 feet/sec 
V.J. = 90 feet/sec 

Substituting these in Equations (55) - (59) yields: 
h, = 17.45 feet 



(63) 


(64) 


Hy = 33,958.57 feet 

k = 18,333.1 feet (65) 

J = 34,758.19 feet . 

Evaluating Equation (62) indicates y is small enough to be assumed zero 
everywhere. The resulting time histories for this example are presented in 
Figure 20. At touchdown the rate of descent is less than 0.2 feet/second 
and the flight path angle is approximately zero. For a more complete 
landing analysis, the investigator might consider the effect of various 
approach and touchdown velocities; the effect of several other glide descent 
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0.0 



Figure 20. Time histories of several parameters 
governing an example landing. 
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angles and possibly, how the vehicle responds for flare at different 
altitudes. These poss ib i I i ties indicate the flexibility with which this 
solution technique may be employed to conduct a landing analysis or a study 
of other performance criteria involving maneuvers. 

Since static predictions are only special cases of the general integrated 
solutions, it is interesting to compare the results obtained by the two 
methods. By specifying variables in a manner similar to the assumptions 
made for derivation of the static equations, the computational accuracy and 
the effects of weight changes can be assessed. 

First, the maximum vehicle speed at any altitude is calculated by simply 
specifying the altitude as a constant and power as the maximum available. 

The integration is terminated once the velocity has settled to a near constant 
value. Results are presented in Figure 21 for various altitudes. 



Figure 21. Determination of maximum velocity as 
a function of altitude. 


At sea level, V max was found to be 254 feet per second which agrees well with 
the static approximations in Table 1. This same result could have also been 
obtained by specifying power as and using a family of curves for differ- 

ent constant velocities. The maximum velocity would be the entry for which 
the altitude remained constant (h = 0). For a velocity greater that V max 
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(level flight) the vehicle would dive and a velocity less than V max would 
induce a cl imb. 

Next, minimum power and the velocity for P m j n were determined at a 
particular altitude. The specified variables were: (1) h equal a constant 

and (2) V as a constant incremented throughout a plausible range. Each 
constant velocity solution produces only one value for P m j n in Figure 22. 


Power 

HR 



Figure 22. Determination of P m j n and speed for P m ; n 


The results are for sea level altitude, but the process is easily repeated 
at other altitudes. As Figure 22 indicates, the minimum power was about 
37.0 hp (20,350 ft-lbs/sec) with a corresponding velocity of 97.0 feet per 
second. These results compare favorably with those of Table 1. 

The maximum climb angle (Ymax^ and +h e velocity for Ymax were also 
calculated. Power was specified as the maximum available and velocity was 
made constant for each of several trajectories. The maximum climb angle 
corresponding to a particular flight path and its specified velocity is then 
plotted on Figure 23. The overall Ymax * s +^en selected to be the largest 
of those in Figure 23. Its associated velocity is then termed as the 
velocity for maximum climb angle. Even though the maximum flight path angle 
(13.15 degrees) occurs at 82 feet per second, any velocity from 65 to 110 
feet per second will produce a Ymax wi+hin one degree of the absolute 
maximum. Thus, the maximum climb angle is somewhat insensitive to velocity 
changes over a rather wide speed range. Even though, this is considered to 
be the maximum climb angle, larger flight path angles may be obtained for 
short periods of time from an oscillatory climb similar to the phugoid 
motion described by Figure 18. However, an average over these oscillations 
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13.20 



Figure 23. Determination of y max and speed for Y max * 

always yields, at most, a mean climb angle with magnitude equal to the Ymax 
from Figure 23. Point performance estimates are similar as indicated in 
Table 1. 

To determine the minimum velocity at an altitude, power was specified 
as the maximum available and velocity was again set equal a different constant 
for each of several paths. The initial altitude is designated as the height 
of interest and the equat ions # are integrated for several minutes to determine 
whether the vehicle ascends (h positive) or descends (h negative). The 
velocity for zero rate of climb is then the minimum velocity for the altitude 
of interest. Figure 24 presents calculation of V m j n at sea level. This 
value of V m j n agrees well with that of Table 1 obtained from point perform- 
ance using a parabolic drag polar. Table 2 illustrates a more realistic 
minimum speed correspond i ng' to an improved drag polar. 



Figure 24. Calculation of minimum velocity for level flight. 
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Finally, calculations were made of absolute and service ceilings. 

Power was specified as the maximum available. A family of 90 minute tra- 
jectories was calculated with each entry corresponding to a different 
specified constant angle of attack. This is not meant to imply that the 
primary objective was to obtain absolute and service ceilings as a function 
of angle of attack. Instead, a was merely a convenient parameter for 
creating a family of curves from which to choose the maximums. Two altitude 
points from each trajectory were recorded in Figure 25. The first was the 
altitude at which h = 1.67 ft/sec and the second was the final height after 
integrating over 90 minutes of real time. The peak in the lower curve of 
Figure 25 is termed the service ceiling and the maximum of the upper curve 
is called the absolute ceiling. Since in actuality the vehicle will continue 
to climb slowly as fuel is consumed, absolute and service ceiling are 
associated with a particular weight. The trajectories represented by 
Figure 25 used approximately 42 pounds of fuel in the 90 minute trajectory. 
The increase of 100 to 150 feet in these values of absolute and service 
ceiling as compared to those of Table 1 is a consequence of the weight 
reduction. The velocities corresponding to the peak points of Figure 25 
are 153.6 ft/sec for the service ceiling and 158.0 ft/sec at absolute 
ceiling. These were likewise increased due to the weight reduction. 

25000 
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h 21000 
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Figure 25. Determination of absolute and service ceilings. 

The preceding discussion of calculation procedures for recovery of 
several static performance parameters indicates that velocity for Ymax is 
not critical and discrepancies in other parameters are related to weight 
variations. 
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The foregoing examples have demonstrated that if desired the path 
performance method can be used to recover the usual static performance 
parameters. In addition, this method can be used to determine the per- 
formance of the aircraft during accelerating flight or in response to a 
complex schedule of pilot control inputs. 


In passing, some comments may be made regarding range performance at 
"practical speeds", i.e. , those obtained at high percents of rated 
horsepower. Table 5 shows that when a < ctL/Dmax ran 9 e ' s reduced. 

In the example given a = 3° rather than 5.85°, a for ( L/D) max . In this 
case, the range is about 5 $ less. This calculation assumes that the 
specific fuel consumption is the same at all power settings and engine 
speeds. As noted in Appendix H the specific fuel consumption may be as 
much as 1.4 times as high at rated power as at the power required for 
maximum range. Thus, not only does the aircraft operate more inefficiently 
but so does the engine. By dropping back to about 75$ power the engine 
can be made to operate very near maximum efficiency so that under these 
conditions one must contend only with aerodynamic inefficiencies. A 
calculation not discussed above was performed for P = P m ax» h = 10,000 ft. 
at the same specific fuel consumption as used in the construction of 
Table 5. The range was only 63/6 of that achieved by flying at a = ot^/Q 
so that at 75$ power one would expect at least a 25$ decrease in n 

maximum range. At maximum power one could reasonably expect a 55$ 
decrease in maximum range. 


max 
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CONCLUDING REMARKS 


The estimation of many light aircraft performance characteristics has, 
as a result of the computer programs provided herewith, been reduced to a 
quick, inexpensive procedure. The programs require that the user supply 
estimates or experimental measurements of the lift and drag variations with 
angle of attack and the variation of power into the airstream as a function 
of speed and altitude. The drag data need not be parabolic. Any accurate 
set of experimental points can be fit precisely by the program and utilized 
in the estimation. Similarly, the program can also fit and utilize arbitrary 
power data. As presently written the program assumes linear variations of 
C|_ with angle of attack and fuel flow with power, but these restrictions are 
not inherent in the method. A knowledgeable user can modify the programs 
without undue difficulty so that they will accept any type of variation he 
might wish to employ. The authors would be pleased to provide suggestions 
on how this might be accomplished to interested users. 

The path performance programs provided with this work offer the user 
an opportunity to study many facets of the flight of light aircraft in 
addition to the equilibrium performance parameters. While the programs do 
not provide a rigorous method of finding the optimum path performance, they 
are sufficiently inexpensive to use that they can usually be employed by 
one with some understanding of the physical situation to arrive at paths 
indistinguishable from the optimal through a trial and error process. It 
is felt that substantial flight time could be saved in evaluating new designs 
by using analytical solutions to indicate the paths of greatest interest. 

Ultimately, one would like to be able to specify just the vehicle 
geometry and its powerplant and have a program which will take this infor- 
mation and compute the vehicle’s performance. By repeating the process 
with some discretion it will be possible to determine a configuration for 
optimum performance subject to the usual constraints of costs, fabrication 
techniques, passenger size, etc. Development of the necessary methods 
and programs for carrying out this procedure is planned for the near future. 
The programs presented in Reference 2 can then be used to compute the 
stability and riding qualities of the optimum performance configuration. 

If necessary, the configuration can be altered and the process repeated 
until the desired characteristics are achieved. 
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APPENDICES 


APPENDIX A - Nomenclature 


AR 
a, a 
BHP 


C 'j' C 2j’ C 3j’ C “J 

C D 


c Dj 

C Do 

C L 


C L(L/D) 


max 


C L(a=0) 


c 


D 


e 

F 

g 

h 

h 


K s 

Kj 

1 o 


L 


» 

L 


aspect ratio 
acceleration 
engine brake horsepower 

coefficients for Spline curve-fit procedure 
drag coefficient 
induced drag coefficient 
zero- lift drag coefficient 
I i ft coefficient 

lift coefficient after landing flare 

lift coefficient for maximum lift to drag ratio 

lift coefficient for zero angle of attack 

specific fuel consumption 

drag or propeller diameter 

Oswald's efficiency factor 

force or take-off weight correction factor 

acceleration due to gravity (32.2 ft/sec 2 ) 

altitude 

aircraft rate of climb 
take-off time coefficient 
take-off ground run coefficient 
static thrust coefficient 
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coefficients of a general drag polar 
lift 

lift after the landing flare 


mass 


N load factor (CjVc^) or engine speed in revolutions 

per minute 

n engine speed in revolutions per second or load factor 

<C»/C > 

P power actually put into the airstream, sometimes called 

thrust horsepower, i.e. , thrust x velocity 

R radius of circular arc traversed in descending from 

an altitude of 50 feet to touchdown 

engine speed in revolutions per minute 

wing area or distance 

thrust or minimum time to climb 

time 

maximum thrust horsepower 
velocity , 

weight 

horizontal distance traveled 
angle of attack 

angle of attack for maximum lift to drag ratio 
sidesl ip angle 
increment in time 

flight path angle in the vertical -horizontal plane 



(see Figure B-1) 

n 

propeller efficiency 

y 

coefficient of friction 

TT 

3.14159 

P 

density of air 

Po 

density of air at sea level 

0 

p/p , ratio of density of air at an altitude to the 
density at sea level 


RPM 

S 

T 

t 

+h Pm 

V 

w 

x 
a 
a 
3 

At 

Y 


(L/D) 


max 
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cj) flight path angle in the vertical -co-normal plane 

(see Figure B-1 ) 

ip flight path angle in the hori zonta I -co-norma I plane 

(see Figure B-1) 

Subscri pts: 


A 

approach 

B 

braki ng 

G 

ground run 

LOF 

lift-off 

max 

maximum 

min 

minimum 

N 

norma 1 to f 1 i ght path 

s 

take-off 

TD 

touchdown 

w 

wind 

zL 

zero lift 

50 

an altitude of 50 feet 

o 

take-off 


A dot over a quantity denotes the time derivative of that quantity. 
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APPENDIX B - Derivation of General Performance Equations 


Those undertaking the estimation of aircraft performance characteristics 
for the first time will usually find it instructive to review the theoretical 
basis of the methodology. If one considers that the word performance 
connotes such things as "how fast will it go?", "how long will it take to 
get to 10,000 feet", etc., he realizes that he is asking questions about 
the motion of the aircraft under certain constraints. The classical means of 
describing the motion of a rigid body in space is through solutions of the 
equations of motion, mathematical statements of Newton's Second Law of Motion. 
It is not surprising that the same equations are also used to study the 
stabi I i tv of the motions, the problem usually termed the dynamics of the 
airplane. The study of aircraft performance, then, can be considered as one 
view of the general problem of aircraft motions while the study of stability 
and control simply views the same problem in another light. The difference 
in the two views, as will become clear from the development below, is 
primarily one of time scale. Stability and control analysis is concerned 
with transient disturbances from an equilibrium motion. The disturbances 
of interest usually have periods of less than 30 seconds. Performance 
analysis on the other hand, is concerned primarily with quasi -equi I ibrium 
flight. Further, the period of interest generally exceeds 30 seconds. 

Thus, the following development proceeds in roughly the same fashion as the 
development of the stability equations. The points of departure and 
difference are noted. 


To begin, make three assumptions: 

Assumption I. The earth and its atmosphere are flat and non- 
rotating. There is no motion of the atmosphere 
with respect to the earth. Accelerations measured 
with respect to a Cartesian coordinate system 
fixed in the earth's atmosphere are therefore 
true accelerations in inertial space. 

Assumption II. The motion of interest is that of the aircraft's 
center of gravity only. Other motions, for 
example rotation about the center of gravity, 
are of interest only insofar as they affect the 
translational motion of the center of gravity. 


Assumption III. While the mass of the aircraft does not remain 
constant with time the change in momentum 
associated with the fuel mass ejected with 
the engine exhaust may be neglected in 
computing the motion of the aircraft. 


As a consequence of these assumptions one may describe the aircraft as a point 
mass and its motion by three equations: 
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ZF X 

= Max 


EF y 

II 

2 

Q) 

'C 

' (B-1) 

ZF 

= Ma 



This is in contrast to stability analysis where the motions of the 
aircraft about its center of gravity are of great interest and require the 
inclusion of the body's three rotational degrees of freedom. On the other 
hand, the time frame treated by stability analysis is sufficiently short that 
the mass can be considered to be constant. 

We take as our coordinate system an x-axis parallel with the earth 
and pointed in the direction of motion at the beginning of the time of 
interest; a y-axis pointed to the right of the direction of motion; and a 
z-axis pointed straight down. In contrast with stability analysis where the 
axis system rotates with the airframe this axis system, once chosen, remains 
fixed with time. In this system the aircraft velocity has three components: 

x = V cos y cos ip 

y = V s i n ip cos y ( B-2) 

h = - V sin y 

Y and \p are flight path angles as defined in Figure B-1 . 

Substition of Equations (B-2) into (B-1) yields 

• • • 

ZF X - M (V cos Y cos ^ - Vy sin y cos ip - Vip cos y sin \p) 

• • • 

EFy = M (V sin ip cos y + Vi p cos y cos \p - Vy sin y sin \p) (B-3) 

• • 

EF Z = M ( — V sin y - Vy cos y) • 

Assumption IV. The forces acting on the aircraft are lift, drag, 
thrust, and weight. The lift acts normal to 
the flight path, the thrust* and drag are 


* Assuming the thrust to act along the flight path is equivalent to 
assuming that the angle of attack is always small. It may therefore seem 
strange that the analysis employs expressions which provide accurate repre- 
sentations of the lift and drag to large angles of attack. This practice is 
justified because (1) light aircraft generally have insufficient power to 
operate continuously at very high angles of attack, (2) light aircraft are 
usually powered by piston engines and unshrouded propellers which makes it 
difficult to assert that the thrust always acts parallel to the wing chord, 
particularly at large angles of attack and (3) it is a simple matter to 
substitute expressions such as T cos a for T and l_' + T* sin a for L in the 
equations without materially affecting the method of solution or the computa- 
tion time. The effect of power at angle of attack is treated in considerable 
detail in section 5.0 of Reference 13. 
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parallel to the flight path, and the weight 

acts along the positive z-axis. There is 

never a side force. This implies that the 

aircraft is symmetrical and is never yawed 

with respect to the flight path (3 is always zero). 

With this assumption, the left hand side of Equations (B-3) can be written 

EF X = (T - D) cos y cos ip - L (sin y cos <j> cos \p + sin ip si n <j>) 

EFy = (T - D) sin ip cos y - L (sin \p sin y cos <p - sin <J> cos ip) (B-4) 

EF Z =(T-D)siny-L cos 4> cos y + |j- . 

if is a flight path angle also defined in Figure B-1. 

Combining (B-3) and (B-4), one obtains 

(T - D) cos y cos ^ “ L (sin y cos <f> cos + sin \p sin <f>) = 

• • • 

M (V cos y c °s ^ - Vy sin y cos ip - ViJj cos y sin i|>) 


(T - D) sin ip cos y - L (sin ip sin y cos <|> - sin <f> cos ip ) = 

. . . (B_5) 

M (V sin ip cos y + V\p cos y cos ip - Vy sin y s i ) 


- (T -D) siny-L cos 4> cos y + W = 

• • 

M ( — V sin y ~ Vy cos y) 

These equations provide the most general description of the motion of 
an aircraft acted upon by L, D, T, and W. Although the aircraft is thought 
of as a point mass, the magnitudes of the lift and drag forces are considered 
to depend upon a, the inclination of the aircraft to the flight path. 

Assumption V. Motions of interest lie entirely within the x-z 
plane. <f) and ^ are therefore zero. 

This of course ignores possible interest in turning flight; however, as a 
consequence of Assumption V, Equations (B-5) simplify easily to 


(T - D) cos y - L siny = M (V cos y - Vy sin y) 

W 

- (T -D) siny- L cos y+W=— (-Vsiny-Vy cos y) 

9 


(B-6) 
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Multiplying the first equation of (B-6) by cos y and the second by -sin Y 
gives 


(T - D) cos 2 
(T - D) sin 2 y + L 


Y - L sin y cos Y = cos 2 y “ V Y sin Y cos y) 
sin Y cos Y “ w s i n Y = g- (^ sin 2 y +Vy cos y sin y^ 


(B-7) 


Adding the two equations one obtains 

(T - D) - W sin y = |v 

Multiplying the first equation of (B-6) by sin y and the second by cos Y 
yields 

? W * • 

(T - D) cos y sin y “ L sin y = — (V cos y sin Y • Vy sin y) 

W 


(B-8) 


(B-9) 


- (T - D) sin Y cos Y “ L cos y + W cos y = " <-V sin y cos Y ~ v Y cos y) 


Adding these two equations, one obtains 

| Vy = L - W cos y- (B- 10 ) 

Equation (B-8) and (B-10) are just simpler equivalents of (B-6). By writing 
the lift and drag forces in the equations in terms of the usual aerodynamic 
quantities and the thrust in terms of the power which is more appropriate 
for piston engine aircraft, (B-8) and (B-10) become finally 

* qP q SV 2 

V = WV • w c D (a) * P ( h) • — - g sin Y (B-11) 

• q SV q 

Y = C L (a) • p(h) *2 y cos y . (B-12) 

The reader will now observe that we have chosen to describe quasi-steady 
aircraft motions by two first-order, non-linear, ordinary differential 
equations. The dependent variables in these equations are P, V, y> W, a > 
and h, while the independent variable is time. Thus to obtain a determinant 
system, one must supply four additional, independent relationships involving 
the dependent variables. One such relationship follows immediately from a 
definition of the rate of climb: 


h = V sin y . (B-13) 

Assumption VI. The fuel flow rate is directly proportional 
to the power developed by the engine. While 
not strictly true, most engines have a region 
around cruise power where the specific fuel 
consumption is nearly constant. The propeller 
efficiency under these conditions is also nearly 
constant. 


74 


(B-14) 


As a result of this assumption one may write a second auxiliary 
relationship: 


W = - cP . 


If necessary, c can be generalized. For turbojet aircraft the fuel flow 
is approximately proportional to thrust output. 

Unfortunately, no other general relationships among the dependent 
variables are known. It is therefore necessary to specify a priori time 
histories of two of the dependent variables in order to obtain unique solu- 
tions to the system of equations. This situation is familiar to pilots who 
realize that the aircraft's trajectory is dependent upon the variation of 
elevator position (speed) and power setting with time. The pilot also 
recognizes that changing the weight and operating altitude also changes the 
power setting and elevator angle one must employ to fly a given path or, 
conversely, for a given power setting and elevator angle, it changes the 
path one flies. With six variables there are 15 different combinations one 
may use to provide the two additional constraints needed. In view of 
Equation (B-14), however, power and weight cannot be specified independently. 
This leaves a total of 14. 

Equations ( B— 11) and ( B— 12) contain in addition three implicit functions 
which must be provided if numerical values are to be obtained. The function 
p(h) is of course the variation of density with geometric altitude. This is 
taken as 

p(h) = p o (1.0 - 6.86 • 1 0~ 6 h ) ** * 2 6 . (B-1 5) 

The functions C^(a) and Cq(cx) depend for their values upon the aircraft 
whose trajectories are desired. Although a third order or higher polynomial 
would be required to represent C|_(a) for alia from - C|_ max to + C Lmax , it 
is usually adequate to employ 

C[_(a) = Ci (a " a zL^ ( B— 16) 

for all speeds above 1.2 Vgj^|_|_p^p^ ^p. Since the system of equations 

must be solved by a forward integration technique in any event, it does not 
add significantly to the overall computational complexity to choose a higher 
order function to represent C|_(a) if this should appear desirable. Similar 
comments can be made with respect toCpta). It is usually represented by 

C D (a) = C Dq + k [C L (a)] 2 (B-1 7) 

although more elaborate descriptions may be used. The addition of a term 
kiCC|_(a)j^2 will permit one to represent the drag coefficient at high angles 
of attack approaching stall quite accurately. 

One final practical constraint must also be observed: the power required 

during any portion of the flight cannot exceed the maximum power available 
for that speed and altitude. Thus if one of the specified variables is not 
power one must always test the computed power required to insure that 
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p £ p, 


( h, V) . 


(B-18) 


max 


Equations ( B— 11) through (B-18), collected below for convenience. 


_ £ 
wv w 


sv" 

2 


Cq(cx) 


p(h) - g sin y 


• g _ . . ... SV g 

Y = w C L (a) • p(h) * — - f cos y 

• 

h = V si n y 

W = - cP 

p(h) = p Q (1.0 - 6.86 • lO^h) 4,26 

C, (a) = Ci (a) + C|_ 

a a=0 

C D (a) = C Dq + k CC|_(a)3 + k-j EIC[_ (ot) II ^ 


(B-19) 


plus the initial conditions of all six dependent variables provide a complete 
description of the motion of the aircraft c.g. in response to a set of forcing 
functions. The forcing functions will take the form of time histories of 
any two dependent variables, for example V and h or P and a. These may be 
chosen arbitrarily, but if one wishes to obtain solutions which represent 
the optimum performance it is necessary that he optimize the form of the 
two variables most appropriate to the parameter being determined. Unfor- 
tunately, it is not possible to prove that the two time histories chosen do 
in fact optimize the particular performance parameter. Only by comparing 
the solutions of (B-19) obtained by applying various physically meaningful 
and realizable time histories of the two most appropriate variables can a 
practical optimum be demonstrated. This solution procedure is discussed 
in detail in the path performance section. 


It may be remarked in passing that while the procedure outlined above 
for determining how to fly the aircraft to obtain the best possible value of 
each performance parameter of interest is the most rigorous known, it is 
practical only with the use of a large, high-speed digital computer. It is 
for this reason that the so-called static or point performance parameters 
came into general use in the years before the availability of electronic 
computers. These procedures, developed from Equations (B-19) assuming that 
the dependent variables do not change with time, are discussed in detail in 
the static performance section. Note, however, that whereas the first two 
equations of (B-19) are differential equations, the static performance 
equations are all algebraic and therefore much easier to solve. Even so, 
evaluation of some of the static performance parameters involves solution 
of a quartic equation or even simultaneous solution of a pair of quartic 
equations. Traditionally, such problems have been solved graphically. In 
the present work the problems are solved numerically through the use of 
digital computer programs. 
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VERTICAL 




APPENDIX C - Point Performance Program 


User Instructions 


The program is written in FORTRAN IV and is designed to run in double 
precision on an IBM 370-165 computer with an average execution time of 1.5 
seconds. This program evaluates static performance of a particular aircraft 
and requires the specification of the following input data: 

(1) The number N of power available versus velocity data points to be 
specified; the control parameter ISUP which designates whether 

the engine is supercharged (ISUP = 1 ) or unsupercharged (ISUP = 0); 
and the reference altitude HREF (feet) which is the altitude at 
which the power versus velocity data points are obtained (for a 
supercharged aircraft HREF must be sea level in this program); 

(2) The N data points of maximum power available PA (foot pounds per 
second) versus velocity V (feet per second) with one data point 
per card, power specified first; 

(3) The four coefficients CDO, CDI1, CDI2, and D of the general drag 
polar which has the form CD = CDO + CDI1*CL 2 + CDI2*CL^, and the 
wing area S upon which the lift and drag coefficients are based; 

(4) The aircraft weight W (pounds), the initial altitude HI (feet) at 
which all performance calculations are made and a final altitude 
HF (feet). If a minimum time to climb schedule and a most 
economical climb schedule from HI to HF, in increments of 100 
feet, are desired then HF must be greater than HI. I f no 
schedule is desired HF must be zero. 

Statements (1) through (4) represent a complete set of data for a particular 
aircraft. Table C-1 gives the format specification for this data. 

Upon completion of the performance calculations with a given set of 
data the program returns to the statement where W, HI, and HF are read. In 
addition to specifying the aircraft weight, the variable W is used as a data 
input control parameter which permits the user to analyze the same aircraft 
for several different values of W, HI, and HF and/or analyze a completely 
different aircraft. The use of W as a control parameter offers the 
fol lowing options: 

(1) If W is positive when the new values of W, HI, and HF are read 
the program yields a new set of performance calculations using 
the original drag polar and power versus velocity curve. For a 
given power curve and drag polar the user may exercise this option 
as many times as desired. 

(2) If W is zero the program returns to the first read statement to 
obtain a complete set of data for a new aircraft. Using this 
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option the user may analyze as many different aircraft as he 
desires. Note that options (1) and (2) may be used together 
(see example discussed below). 

(3) If W is negative the program terminates. Thus, the final data 
card for any computer run must have a negative value for W. 

The example data set given in Table C-2 utilizes all three of the options 
presented above. Using this data the program calculates two sets of per- 
formance characteristics for each of the two aircraft specified. A sample 
output for the Cessna 182 is presented after the program listing. 



Table C-1. Format specification for input data. 
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HI* 0.000 


HF* 0.000 


LAST CARD 


it m 1 » i i imiiiiiniimiiiiiiiinniBiiniMiiiiMDmumii c4>44««47aaasiuu$iaMS7M*atie(iNaatiBa*7i7; 


W* 2750.000 


HI* 0.0D0 


HF* 0.000 


i i i 4 s i ) i t»)imiHttw»»ii»naowa»n»nMHBPM»»g«»««!c<m « « i? a wmi ns m ebb a a mi snut,mvmmnnnnt*n*»n 


W* 2500. 000 


HI* 0.000 


HF* 0.000 


$* 180.000 


i i m i i i i iioni 7 »Hnn»B»aaaB*aa*?a»aii»m 4 B;n 7 BaB 4 i 4 i 


COO* 0.04700 CDI1* 0.0ID0 CDI2* 0.044D0 0= 5.4D0 


i t M i i 7 1 1 n u it iik is a it it nan sanaati a nan sum Bapaaati «4i a aait a aasi uu annua 


Bummiounnunnnnnx 


f (PA) 12 - 125000. 0D0 V 12 * 375. 0D0 

il 14 s i 7 i iioiiiiiiKiiiiinminainiiiiiiiiiiiinuMBmiiBnimnimHMMMiMummuiMinummiiiiiMiiimni 


i t 14 i i ? i imiuiiMi)iniliB»imi»B»n»7i»iiiniiiiiiiii»i4nuiiMimtng»n»»»«»iii emtfmum»nnTin*nmrnnm 


(PA) 2 - 40000.000 V 2 = 85.000 


CPA) ^ 0.000 


V^* 0.000 


i t il i i i * >wnum 4 Bwi 7 i«>ianaawaanwnaiiBMa»»t 7 aaa 4 i« 4 i 44 «a 47 aa»uuu> 4 »awa»>n«tMwaaya)on 7 tnM»ai 7 m 


H» 12 ISUP- 1 HREF* 0.000 


i t ms i 7 • »Miiitin 4 BMi 7 B»aflBaMBanaaM]ij;i)MBMnaBa 4 iaua«« 4 r««a»»UMMMsiii 


Hunan nnnMnit >7 naio 


W* 0.000 


HI- 0.0D0 


HF* 0.000 


i i il i * 7 i > a ii » nw is w if a a n a a na a r? a n » n w a agjii? a a a n 4t u a a a <7 a a a si « a a Hawaii a a a a a a a it m m a n nTiHTsannna 


W- 2900.000 


HI- 0.000 


HF* 5000.000 


* i i ? i iwiiuut4Bai?aiiatiaawaat7ana«aaMa»pa 


a sia a to nnnnnn}innn 


W- 2650.000 


HI* O.ODO 


HF* 5000.000 


i i ii i > 7 i ianitiiii5ai7aiiatiaawaai;aBana8waawaaaiicua«ai?a4tasiMusiaai7BaaiiwaMaai7aawii77Brtaa7?an» 


S- 174.000 


r;ana]iBas»»wi 7 aaaiiawwaa 47 aaa>iuu»a»s 7 M»auBuwaas 7 aaan 7 t 7 i/ 4 Ba 7 tana 


COO- 0.026900 C0I1* 0.05400 CDI2* 0.01800 D- 6.500 


i taiittu»aaii»BaBBa»a»rM 


a 4 i«aMaa 47 aaaii»UMaat 7 aaa«BaHaawaa»nnn/ 4 aniin 7 *B 


(PM 15 - 100000.000 V,,* 400.000 

i t i i l l i i imniiiHHiiiri nan an a aa v a a a n BBMaaMaaanaaaaaitaaawHuwawwaaaaim 


it I 4 ii ii ia!iituK»ai 7 aaanaaNBat 7 aaaanaNBBiiaaa 4 i* 4 iMnMiiaaaM»u* 4 Biis 7 Baa«en«aai 7 aa* 7 i 7 ! 7 ) 74 niir 7 a 7 ia 


(PM 2* 50000.000 V 2 * 80.0D0 


i i i i l * t i taimunnnitaaanBnwaanaaBiiBaaaBgaaaiiituaaatiaaasiHUMawstaaatiBUMaatiaaniiTtTtMTinitnna 


(PA) j* 0.000 V,* 0.000 

i 1141 i i i i aiittUHBWuawanaaaaanaaan aaManttnaaiiauMnauanasistusiBau aaasi CO want) a a »?i it ti * is» » ana 


N* 15 ISUP* 0 HREF* 0.000 


FIRST CARD 


i i t ill i i taHt 20 H»Muaaati&a 


• Mill II I III 1 1 1 II Hill 

n t n m n ii 1 1 1 m m i 


nmimnimmmt 

mimnunminm 


I Mill 1 1 4 4 4 1 4 4 4 4 4 4 4 4 4 4 
IIHIHHHHHHH5M 


1111111111111111111111 

miinmiiiininin 


1111111111111111111111 

iiiiiiiiiiitmiiiiiii 



1 1 1 1 1 1 • 1 1 1 1 1 1 1 1 1 1 1 1 1 1 • 1 1 1 1 1 1 • 



1 1 1 1 FVKitMr 7 mm 

COMPUTING CENTER 


llllllllll IIIIIMM 

MMMMMMMMMMMI MMM 


Table 0,-2. Example data set for point performance program. 


N - NUMBER OF POINTS ON PA VS. V CURVE 
ISUP - CONTROL VARIABLE FOR SUPERCHARGEO ENGINE 
ISUP - 0 -> ENGINE NOT SUPERCHARGED 


•c-o®'0««®'«r-f*»h-p»r'-r-f.f>-r»r'®®««>®*>*)®«®0'ff‘ff'0'a»o-0'«r«roooocoo©eo — — »-■ — — ****** nnnnnn 


o> 

c 



«M 


fin<NB9'0-HN'n4«'4N®e'0 
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TION POLY FCR DEFINITION OF PSUEOO-POi YNON( HI 14 C 

WHOSE ROOTS ARE VMAX ANC VNIN 15 C ONLY ONE ROOT HAS BEEN FOUND. DETERMINE IF IT IS VNIN OR VHAX 
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VRNGHX - VELOCITY FOR HAXINUH RANGE 

COEF - COEFFICIENTS OF RSUE DO-POLY NON I AL WHOSE 
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CLIMB AND COMPARE IT TO RATE OF CLIMB OF ANY PREVIOUS ROOT 
TO DETERMINE MHICM IS THE GREATER. 
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IF VV IS LESS THAN Vtl». EXTRAPOLATE TO VV USING PACOEF’S 26 C 

THAT ARE FOR VELOCITY RANGE V(l» TO V!2». 27 FUNCTION POLV(N,A,XI 

28 REALMS A ( 1 3 ) • POL Y • POL V 2 * X 
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35 C N NUMBER OF TABULATED DATA POINTS 

FI « POL Y (NOIN* COEF » X 1 ) 36 C Y TABULATED FLECTION VALUES 

IFIFI.NE.O.OOOt GO TO 1 37 C X TABULATED ARGUENENT VALUES 

XRCOT - XI 3B C AA ARRAY OF COEFFICIENTS OF THE CUBIC POLYNOHIAL 
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Sample Output 


POWER AVAILABLE VS* VELOCITY 
REFERENCE ALTITUDE * 0.0 FEET 

PAI FT-LBS/SEC ) VI FT/SEC I 


0.0 


0.0 


0.291500 

05 

0.273300 

02 

0 .524700 

05 

0.546700 

02 

0.699600 

05 

0.820000 

02 

0.8 1 6200 

05 

0.109330 

03 

0.8T450D 

05 

0.136670 

03 

0.90948D 

05 

0.164000 

03 

0.944460 

05 

0.19133D 

03 

0.95B65D 

05 

C.21867D 

03 

C.96778D 

05 

0.246000 

03 

0.973610 

05 

0.273330 

03 

0.97944D 

05 

0.300670 

03 

0.9947CD 

05 

0.328000 

03 

0.9947CD 

05 

0.355330 

03 

0.99470D 

05 

0.382660 

03 


AIRCRAFT CHARACTERISTICS 

CO « 0 • 268800-0 1 • 0.542420-01*CLP*2 ♦ 0. 17751D-01*CL** 0.650000 01 
WING AREA - 0.174C0D 03 SQ.FT WEIGHT - 0.265000 04 LBS 


STATIC PERFORMANCE AT AN ALTITUDE - 0.0 FT 

WITH MINIMUM TIME AND MOST ECONOMICAL CLIMB SCHEDULES TO A FINAL ALTITUDE - 0.100000 05 FT 

++++++++++ ++++* 4 +++++++++++++++ +++++++++++++++++ 4 + 4 + +++++++++++++++++++++++++ 4 + 4 + +++++++++ 

MINIMUM LEVEL FLIGHT SPEED - 0.904650 02 FT/SEC 
LIFT COEFFICIENT - C. 156380 01 DRAG COEFFICIENT - 0.484180 00 

MAXIMUM LEVEL FLIGHT SPEED - 0*252570 03 FT/SEC 
LIFT COEFFICIENT - 0. 200620 00 DRAG COEFFICIENT * 0.290630-01 


MAXIMUM CLIMB ANGLE - 0.10220D 02 DEG 
VELOCITY FOR MAXIMUM CLIMB ANGLE - 0.11710C 03 FT/SEC 
LIFT COEFFICIENT - C. 933300 00 ORAG COEFFICIENT - 0.854750-01 


VELOCITY FOR MAXIMUM ENOURANCE - 0.125710 03 FT/SEC 
POWER FOR MAXIMUM ENDURANCE - 0.275450 05 FT-LBS/SEC 
LIFT COEFFICIENT - 0.809850 00 DRAG COEFFICIENT - 0. 669610-01 


VELOCITY FOR CLASSICAL MAXIMUM 
LIFT COEFFICIENT - C. 634230 00 


RANGE - 0.14205D 03 FT/SEC 
DRAG COEFFICIENT - 0. 496190-01 


SERVICE CEILING - 0.19442D 05 FT 
VELOCITY AT SERVICE CEILING - 0.17554D 03 FT/SEC 
LIFT COEFFICIENT - 0.76424D 00 DRAG COEFFICIENT - n. 616530-01 

ABSOLUTE CEILING - 0.212360 05 FT 
VELOCITY AT ABSOLUTE CEILING - 0.180230 03 FT/SEC 
LIFT COEFFICIENT - C.T7C56D 00 DRAG COEFFICIENT - 0.62348D-01 


MAXIMUM RATE OF CLIMB SCHEDULE FROM 0.0 FT TC 0.100000 05 FT 


HIFTI 


R/C (FT/SEC 1 

VIFT/SEC 1 

PIFT-LBS/SEC* 

CL 


CO 

TISEO 


0.0 


0.222570 

02 

0.136010 

03 

0.873540 

05 

0.691830 

00 

0.544600-01 

0.0 


0.500000 

03 

0.21659D 

02 

0.13662D 

03 

0.85921D 

05 

0.695810 

00 

0.548210-01 

0.227750 

02 

0.100000 

04 

0.21065D 

02 

0.13725D 

03 

0.845040 

05 

0.69966D 

00 

0.55 1 74D-01 

0.461850 

02 

0.150000 

04 

0.20475D 

02 

0.137090 

03 

0.831020 

05 

0.703390 

00 

0.555190-01 

0.70263D 

02 

0.200000 

04 

0.198890 

02 

0.138560 

03 

0.817160 

05 

0.70700D 

00 

0.550560-01 

0.950420 

02 

0.250000 

04 

0.19307D 

02 

0.13925D 

03 

0. 80 3450 

05 

0.710490 

00 

0 .561 85D-0 1 

0.120560 

03 

0.300000 

04 

0.187280 

02 

0.139960 

03 

0.789890 

05 

0.713840 

00 

0.565040-01 

0.14686D 

03 

0.350000 

04 

0.181540 

02 

0 .140690 

03 

0.776470 

05 

0.717080 

00 

0.568140-01 

0.17398D 

03 

0.40000D 

04 

0.175830 

02 

0.14145D 

03 

0.763200 

05 

0.720180 

00 

0.571150-01 

0.201970 

03 

0.450000 

04 

0.170150 

02 

0.142220 

03 

0.75C08D 

05 

0.723160 

00 

0.574050-01 

0.2308 BD 

03 

0.500000 

04 

0. 1645 ID 

02 

0.143020 

03 

0.73709D 

05 

0.726000 

00 

0.57684D-01 

0.260770 

03 

0.550000 

04 

0.158910 

02 

0.143840 

03 

0.724250 

05 

0.728720 

00 

0.579530-01 

0.291700 

03 

0.60000D 

04 

0.153350 

02 

0.14468D 

03 

0.711550 

05 

0.731310 

00 

0.582110-01 

0.323730 

C3 

0.65000D 

04 

0.14782D 

02 

0.145540 

C3 

0.69899D 

05 

0.733760 

00 

0.58457D-01 

0.356950 

C 3 

0.700000 

04 

0.142330 

02 

0.146430 

03 

0. 686S7D 

05 

0.736080 

00 

0.586910-01 

0.391420 

03 

0.750000 

04 

0.136870 

02 

0.14734D 

03 

0.67428D 

05 

0.738280 

00 

0.589 14D-0 1 

0*427260 

03 

0.800000 

04 

0.131440 

02 

0.148280 

03 

0.662130 

05 

0.740340 

00 

0.591240-01 

0.464540 

03 

0.89000D 

04 

0.126060 

02 

0.149240 

03 

0.6501 ID 

09 

0.742270 

00 

0.593220-01 

0.503390 

03 

0.900000 

04 

0.120700 

02 

0.15022D 

03 

0.638230 

05 

0.744070 

00 

0.5950BD-01 

0.543940 

03 

0.950000 

04 

0.115390 

02 

0.151220 

03 

0.626490 

05 

0.74573D 

00 

0.596810-01 

0.586320 

03 

0.100000 

05 

0.110100 

02 

0.152260 

03 

0.614880 

05 

0.74727D 

00 

0.598410-01 

0.630690 

03 




MOST 

ECONOMICAL RATE 

OF CLIMB SCHEDULE 

FROM 0.0 


FT TO 0.100000 05 FT 


HIFT) 


R/C I FT / SEC I 

V I FT /SEC ) 

PI F T-L8S/SEC 1 

CL 


CO 

0.0 


0.221090 

02 

0.128970 

03 

0.862280 

05 

0.769400 

00 

0.6 2 2 200-01 

0.500000 

03 

0.215140 

02 

0.124860 

03 

C. 848770 

05 

0.770160 

00 

0.623040-01 

0. 100000 

04 

0.209220 

02 

0.130750 

03 

0.835360 

05 

0.770890 

00 

0.623850-01 

0.150000 

04 

0.203320 

02 

0.13166D 

03 

0.62205D 

05 

0.771600 

00 

0.624640-01 

0.200000 

04 

0.147450 

02 

0.132580 

03 

0.808650 

05 

0.77227C 

00 

0.625390-01 

0.250000 

04 

0 * 191600 

02 

0.133510 

03 

0.795750 

05 

0.772910 

00 

0.626100-01 

0.300000 

04 

0.185770 

02 

0.134450 

03 

0.782750 

05 

0.77352D 

00 

0. 626780-01 

0.350000 

04 

0.179970 

02 

0.1 354 ID 

03 

0.769860 

05 

0.774090 

00 

0.627430-01 

0.400000 

04 

0.174200 

02 

0.136 390 

03 

0.757080 

05 

0.77463C 

00 

0.628030-01 

0.450000 

04 

0.168450 

02 

0.137370 

03 

0.744400 

05 

0.775120 

00 

0.628580-01 

0.500000 

04 

0.162720 

02 

0.138370 

01 

0.731840 

05 

0.775590 

00 

0.629110-01 

0.550000 

04 

C. 157020 

02 

0.139380 

03 

0.719380 

05 

0.776030 

00 

0.629600-01 

0.600000 

04 

0. 15135C 

02 

0.140410 

03 

0.707030 

05 

0.776430 

00 

0.630050-01 

0.650000 

04 

0.145700 

02 

0.141450 

03 

0.694600 

05 

0.776790 

00 

0.630470-01 

0.700000 

04 

C.140C70 

02 

0.142510 

03 

0.662670 

05 

0.777120 

00 

0.630840-01 

0.750000 

04 

0.134470 

02 

0.143590 

03 

0. 670660 

05 

0.777410 

00 

0.631160-01 

0.800000 

04 

0.126900 

02 

0.14468D 

03 

0.658760 

05 

0.777650 

00 

0.631440-01 

0.850000 

04 

C. 123350 

02 

0. 145780 

03 

0.646970 

05 

0.777850 

00 

0.631660-01 

0.900000 

04 

0.117820 

02 

0.146910 

03 

0.635300 

05 

0.777990 

00 

0.631830-01 

0.950000 

04 

0.112320 

02 

0 .148050 

03 

0.623750 

05 

0.778090 

00 

0.631930-01 

0.10000D 

05 

0.106840 

02 

0.149210 

03 

0.612310 

05 

0.778130 

00 

0.631980-01 


maximum r/c, power available, 
AT 0.0 

R/CIFT/SEC) P A I FT-LB 5/SEC I 


0.391050-09 

0.74225D 

05 

0.136260 

02 

0.7834 30 

05 

0.191130 

02 

0.81825D 

05 

0.213400 

02 

0.B4449D 

05 

0.22151D 

02 

0.864040 

05 

0.222170 

02 

0.879200 

05 

0.218260 

02 

0.892030 

05 

0.211210 

02 

0.904320 

05 

0.201930 

02 

0.917670 

05 

0.190230 

02 

0.931410 

05 

0.175210 

02 

0.943170 

05 

0.156070 

02 

0 .951000 

05 

0.133030 

02 

0.955690 

05 

0.106710 

02 

0.959090 

05 

C. 774840 

01 

0.562540 

05 

0.452220 

01 

0.965920 

05 

0.96874C 

00 

0 « 968910 

05 

0.991510-10 

0.969590 

05 


c POWER REQUI RED VS VELOCITY 
FT 

PRQIFT-LBS/SECI VIFT/SECI 


0.742250 

OS 

0.904650 

02 

0.422340 

05 

O.IOOOOC 

03 

0.311760 

05 

0.110000 

03 

0.278990 

05 

0.120000 

03 

0.277040 

05 

0.130000 

03 

0.290440 

05 

0.14000C 

03 

0.313650 

05 

0.150000 

03 

0.344610 

05 

0.160000 

03 

0.382560 

05 

0.170000 

03 

0.427300 

05 

0.180000 

03 

0.478870 

05 

0.190000 

03 

0.537410 

05 

0.200000 

03 

0. 603150 

05 

0.210000 

03 

0.676320 

05 

0.220000 

03 

0.757200 

05 

0.2 30000 

03 

0.846080 

05 

0.240000 

03 

0.943240 

05 

0.250000 

03 

0.969590 

05 

0.252570 

03 


APPENDIX D - Path Performance Program 


User Instructions 


The program is written in FORTRAN IV and is designed for execution in 
double precision on an IBM 370/165. For convenient description, the program 
is divided into five parts as follows: 

(1) Main I ine - This section handles overall program control, roads all 
input data, converts the variable units so that they are consistent 
for execution but convenient for input and output. Based on the 
initial input, it adjusts the value of sea level density so that 
for the first point lift equals weight to within maximum machine 
accuracy. This is necessary for subsequent calculations to yield 
accurate values of y and y. If an adjustment in p Q of more than 

5 % is needed, the input data is considered inconsistent and 
execution stops. Therefore, the user should always adjust the 
input data so that lift and weight are initially equal. The 
mainline also prints out the integrated solution and tests the 
results for variables exceeding upper or lower limits. 

(2) Subroutine SPLINE - This subroutine provides maximum power available, 
as a function of a h and V, for use in testing the calculated 

power or possibly as an input when power is specified as P max . A 
detailed description of the spline procedure and calculation of 
P m a X (h>V) is given in the section entitled Computerization 
Procedure for Point Performance. 

(3) Subroutine F - Corresponding to a particular pair of specified 
variables, this subroutine calculates both derivatives of the 
variables to be integrated and updates the algebraic parameters 
at each integration point. 

(4) Thirteen FUNCTION Subprograms - These subprograms supply values 
for the specified variables and their derivatives throughout the 
integration region. The specified variables may be functions of 
one or more different flight parameters. 

(5) Subroutine TRENOR - This subroutine integrates the general equations 
using a modified Runge-Kutta Pred i ctor-Corrector technique which 

is described both in the Path Performance section and in Appendix 
G. This subroutine also adjusts the integration step size by 
halving or doubling on the basis of an error criterion. 

Throughout this program, the subscripted variable Y is used to denote 
any of the five integrated variables. Their positions are defined as 
fol lows: 


94 



Y ( 1 ) = range, x 

Y(2) = weight, W 

Y(3) = altitude, h 

Y(4) = flight path angle, y 

Y(5) = velocity, V 

Once a suitable set of specified variables has been chosen, the user 
selects the associated Key number from Table D-1 . 


Key 

Specified Variables 

FUNCTION Subprograms Utilized 

1 

h, V 

H, DH, DDH, V, DV 

2 

h, Y 

H, DH, DDH, GAM, DGAM 

3 

h, a 

H, DH, DDH, ALPHA, DALPHA 

4 

h, W 

H, DH, DDH, W, DW 

5 

h, P 

H, DH, DDH, P 

6 

V, y 

V, DV, GAM, DGAM 

7 

V, a 

V, DV, ALPHA 

8 

V, W 

V, DV, W, DW 

9 

V, P 

V, DV, P 

10 

Y, a 

GAM, DGAM, DDGAM, ALPHA, DALPHA 

11 

Y, W 

GAM, DGAM, W, DW 

12 

Y» P 

GAM, DGAM, P 

13 

a, W 

ALPHA, W, DW 

14 

a, P 

ALPHA, P 


Table D-1. Relation between Key numbers, specified variables, 
and FUNCTION subprograms. 


Based on the Key numbers this program selects the variables which will be 
integrated and those to be specified. Corresponding to a Key number, the 
user must provide relations for the associated FUNCTION subprograms as shown 
in Table D-1. Even though each specified variable and its derivatives are 
implicitly functions of time (the independent variable), they may be explicit 
functions of other parameters as well. If this is the case, it is mandatory 
that the argument lists of both the FUNCTION subprograms and of the calling 
statement in subroutine F be made compatible. The following two examples 
will i I lustrate this procedure. 

Example (1). Choose velocity and angle of attack as the specified 

variables. From Table D-1, this pair has a Key value of 
seven. For Key = 7, Table D-1 indicates utilization of 
subprograms V, DV, and ALPHA. Let velocity equal 151 
feet per second and angle of attack equal that for maximum 
lift to drag ratio, as in Case (3) of Table 5. The 
FUNCTION subprograms now appear as: 
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FUNCTION V(T) 

IMPLICIT REAL*8(A-H, O-Z) 

V = 151. DO 

RETURN 

END 

FUNCTION DV(T) 

IMPLICIT REAL*8(A-H,0-Z) 

DV = O.ODO 

RETURN 

END 

FUNCTION ALPHA (T) 

IMPLICIT REAL#8 ( A-H , O-Z ) 

ALPHA = . 10258D0 

RETURN 

END 

The FUNCTION calls of subroutine F in the Key = 7 section 
would then read: 

29 Y(5) = V(T) 

ALPHA 1 = ALPHA (T) 

PI = Y ( 2 ) *Y ( 5 ) #DV ( T ) /G+ • • • 

Example (2). This example is more complex and is designed to help 

demonstrate the program’s flexibility. Let velocity and 
flight path angle be the specified variables. Table D-1 
indicates Key = 6 and the needed FUNCTION subprograms are 
V, DV, GAM, DGAM. Let V, V, y, and y be those for the 
landing analysis described by Equations (54)-(59). The 
variables are as follows: 



-.04363 radians 


- tan -1 


( X - 34758.19 ) 
H - 18333.1 


17.45 < H < 1500 


0 < H < 17.45 


V 



for all H 
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140.0 


17.45 < H < 1500 


) H + 90. 0 < H < 17.45 

17.45 
and since 


dV = dV dh 
dt dh dt 


50 • = 50 

17.45 h 17.45 


sin y 


Thus, 



0.0 


50 

17.45 


V 


si n y 


17.45 < H £ 1500 


0 < H < 17.45 . 


The FUNCTION subprograms corresponding to each of these 
variables are as follows: 


FUNCTION GAM(X,H) 


IF(H. LE. 1 7.45D0) GO TO 1 
GAM = -.0436300 
RETURN 

1 GAM = -DATAN((X-34758.19D0)/(H-1 8333. IDO)) 
RETURN 
END 

FUNCTION DGAM(T) 


DGAM = O.ODO 

RETURN 

END 

FUNCTION V(H) 


IFOH.LE. 17.45D0) GO TO 1 

V = 140. DO 

RETURN 

1 V = 50 . D0*H/ 1 7 . 45DO+90 . DO 
RETURN 
END 
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FUNCTION DV (H, V ,GAM) 


IF(H.LE. 17.45DO) GO TO 1 

DV = O.ODO 

RETURN 

1 DV = 50 . DO*V*DS I N (GAM) / 1 7 . 45D0 
RETURN 
END 


For these variables the FUNCTION calls of subroutine F 
in Key = 6 section should read as follows: 

28 Y(5) = V(Y(5) ) 

Y(4) = GAM ( Y ( 1 ) , Y ( 3 ) ) 

CL = (DGAM(T)+G • • . 

PI = Y(2)*Y(5)*DV(Y(3),Y(5),Y(4))/Gt • • • 

These two examples in addition to those of the sample program, which appears 
later in this Appendix, serve to illustrate the steps necessary for supplying 
the specified variables. 

For analyzing a particular aircraft, this program requires input data 
which may be grouped into three categories. Each input variable is defined 
in detail on the first page of the program listing. 

Group (1). The first group, composed of all information necessary to 
calculate maximum power available as a function of velocity and altitude, is 
identical with the initial set of data cards for the point performance 
program. The number N of power available versus velocity data points to be 
specified; the control parameter ISUP which designates whether the engine is 
supercharged (ISUP = 1 ) or unsupercharged (ISUP = 0) and the reference 
altitude HREF (feet) which is the altitude at which the power versus velocity 
data points are obtained (for a supercharged aircraft HREF must be sea 
level in this program) are read from the first data card. The next N cards 
contain data points of maximum power available PA (ft-lbs per sec) versus 
velocity VA (ft per sec) with one data point per card, and power specified 
first. Thus the first N + 1 cards of the two programs have an identical 
purpose. 

Group (2). The second group consists of variables which control the 
integration of the path performance equations. The value of KEY indicates 
those variables which will be integrated and those to be specified. PRINT 
gives the frequency of print out. MAXHLV limits the net number of times for 
halving the step size. ACCHLV and ACCDBL govern the halving and doubling 
of the step size and HMAX is the maximum step size permitted. 

Group (3). The third group includes initial conditions for the 
trajectory, pertinent aircraft character ist ics, the maximum time for the 
trajectory, and other necessary parameters. These variables are defined, 
complete with units, in the program listing. 
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Table D-2 gives the format specification for the input of all the 
necessary data. 



Table D-2. Input format specification for data of the 
Path Performance Program. 


This program is designed so that the data of Group (1) (for computing 
p max (h * v ^ is read on ly once during the program execution. For a particular 
vehicle with P max described by Group (1) data, sets of data from Groups (2) 
and (3) may be repeatedly read for calculating several different trajectories. 
Upon completion of a single trajectory, the program returns to the statement 
where KEY, IPRINT, . . . are read. If KEY is any integer between 1 and 14 
the program reads a new set of Group (2) and (3) data. This process may be 
often repeated. After the last trajectory a card with KEY equal zero is 
inserted to stop the program execution. 


The 

appendix 


program listing and sample output which appear at the end of this 
are used to illustrate the integration of two different trajectories 
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during one program execution. These two examples will adequately illustrate 
the program's operation, even though several other cases could have easily 
been included. Compatibility of the FUNCTION subprogram argument lists with 
their calling statements may be observed in the program listing. The sample 
output resulted from the program's execution of the example data set in 
Table D-3. The following discussion describes collection of input data for 
each example with the total data set presented in Table D-3. The airframe 
for these examples is presented in Figure 2, and the power (PA) versus 
velocity (VA) points were obtained from Appendix F. 

Example (1). An aircraft initially weighing 2650 pounds and having a 
parabolic drag polar is to fly with angle of attack for 
maximum lift to drag ratio at a constant altitude of 
10000 feet. 


Let 

c L 2 

C D " C D 0 + ireAR 


.0269 + 


7r( .98) (7.378) 


thus 

ot( l/d ) max = • 10258 radians = 5.8778 degrees. 
Also, 


C, = Ci + C, a = .309 + 4.608( . 10258) = .7817 

L (a=0) L <* 

and 


C D = .0538 . 

At the initial point lift must equal weight, so 
±p(h 0 )V*S C L = W Q 
or 


2W 


Vo = ( 


p(h 0 )C L S 


) 2 = 


2(2650) 


( . 00 1755) ( .7817) ( 174) 


i_ 

•) 2 


V 0 = 149.0 ft/sec 
and 


p o = ±p(h o )C 0 v£ = 27173.0 ft-lbs/sec = 49.4 hp . 
Thus the initial conditions are as follows: 
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h Q = 10000. 
V 0 = 149. 

Yo = 0.0 
a Q = 5.8778 


W Q = 2650. 
P 0 = 49.4 
x Q = 0.0 


Values for the other necessary parameters are as 
fol lows: 


c = 0.6 
TDEL =0.5 
TMAX = 480.0 
G = 32.2 
S = 174. 
p Q = .00238 

C, = 4.608 
L a 


C Do = .0269 
E = .98 
AR = 7.378 
VMIN =35.0 
WEMPTY = 2150. 
EX = 2.0 

C L(a=0) = 0,309 


For a and h specified. Table D-1 indicates that KEY = 3. 
The correct FUNCTION subprograms associated with this 
value are shown in the program listing. Typical values 
of the other parameters are as fol lows: 

IPRINT = 10 ACCHLV = .001 

MAXHLV = 2 ACCDBL = .00002 


Since these inputs produce rather smooth solutions a 
maximum step size of 60 seconds (HMAX =60.) is employed. 

The second example will consider an airframe slightly different from 
that of the first example but described by the same maximum power curves. 

Example (2). A 2700 pound aircraft initially flying with velocity of 
120 ft/sec at sea level and having a non-parabol ic drag 
polar is to climb under full throttle for 30 minutes with 
a constant flight path angie of 1.5 degrees (.026178 
radians). Let 

C D = .02987 + .07276 Cl*" 1 

then 

C Do = .02987 
and 

• TTeAR = - 07276 or e = 17(7. 378) (.07276) ,593 * 


Since lift and weight must be made equal at the first point, 
the initial angle of attack is computed as follows: 
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2W 0 2(2700) 

C L = p(h 0 )V 0 2 S “ (.00238) (120)* 174 = - 9055 

and 

a 0 = (Ci - C L/ _.)/C L = (.9055 - .309)/4.608 

a 0 = .12945 radians = 7.42 degrees 


P 0 = p max(h 0 ,V 0 ) = P max (0.0, 120) 
The complete initial conditions 


VQ “ I • 

Yo = 1.5 
a 0 = 7.42 

Other necessary 

c = 0.6 

TDEL = 0.5 

TMAX =30.0 

G = 32.2 

S = 174. 

p Q = .00238 

C, = 4.608 
La 


parameters are: 

c Do = 
E = 

AR = 
VMIN = 
WEMPTY = 
EX = 
C, 

L (a=0) 


= 153.5 hp 

are as to I lows: 

2700. 

153.5 

0.0 


.02987 

.593 

7.378 

35.0 

2150. 

2.991 

0.309 


For y and P as specified variables. Table D-1 indicates 
KEY = 12. The correct FUNCTION subprograms associated 
with this value are shown in the program listing. Typical 
values of the other control parameters are as follows: 

IPRINT = 5 ACCHLV = .001 

MAXHLV = 2 ACCDBL = .00002 

Since flight path angle is restricted to be a constant, 
no oscil I at ions are present and HMAX is set equal to 30 
seconds. The input data for these two examples is 
presented in Table D-3. The solution time histories for 
these two examples are presented in the sample output. 

The overall program has some limitations which should be brought to 
the user's attention. Several of these could be overcome by increasing the 
program’s complexity, but for the analysis of light aircraft this was deemed 
unnecessary. The lift coefficient, C|_, was restricted to lie between 0 and 
15 so as to prevent the inclusion of a general root solver. A maximum flight 
path angle of one radian was imposed solely as a check point. For most Key 
numbers the program may be used for flight path angles having magnitudes of 
nearly 90 degrees. In the derivation of the general equations of motion, 
power was assumed independent of angle of attack. If this restriction is 
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LAST CARD 

I ? 14 5 i 7 t 9 10 11 12 1 J It IS K I? U 19 3 21 72 3 W 25 X 27 3 3 30 11 12 11 M * X 17 1 H 4 # 41 47 41 M 4 ^( 7 (ll 9 M 5 l 52 53 i (55 54 5 7 5 l 59 miKilHia<H 7 Mai 0 7 l 77 7 3 1 ( 75 ?t )7 7 S) 9 «) 


' E 

AR 

VMIN 

WEMPTY 

EX 
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3.5D+01 
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TMAX 

G 

s 
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CLA0 

CLA 

CD0 ^ 
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3.D+01 

3.22D+01 
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2.38D-03 

3.090-01 

4.608D0 

2. 987D-02 
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f HI 

VI 

GAM1 

ALPHA1 

W1 

PI 

XI 

> 

C 

O.ODO 

171151719 10 1 

1.2D+02 

1 17 11 K 15 II 17 II 19 70 71 

1.5D0 

27717(7571 77 3310 

7.42D0 2.7D+03 
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1.535D+02 O.ODO 
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0.600 
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f KEY 

IPRINT 

MAXHLV 

ACCHLV 

ACCDBL 

HMAX 


> 

12 

17 1 « 5 1 7 1 9 10 1 

5 

1 1? 11 K 15 II 17 II 19 70 71 

2 

77 7)»I5 7i71»71I 

1.D-03 

3H7U31fiX17X9« 

2.D-05 3.D+01 

lieo44et*(7 3«jlSIBS154BS4S7SIS4«il42U*4ii«4?B3l0 71 72 »/17S?4/7 7|74l0 

f E 

AR 

VMIN 

WEMPTV 

EX 



> 

0.9800 

1 7 1 t 5 1 7 1 9 10 1 
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2. DO 

4l4?4jM4S«l7a«»5lH53!XB44 57B»Bllfi434iaB67Sai0 7l 77 71 <4 15X 77 7l3» 

f TDEL 

TMAX 

G 

s 

RHO 

CLA0 

CLA 

'CD0 

0.5D0 

4.8D+02 

3.22D+01 

1.740+02 

2.38D-03 

3.09D-01 

4.608DO 

2. 69D-02 
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f HI 
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ALPHA 1 

W1 

PI 

XI 

C 

1.D+04 

1.49D+02 

O.ODO 

5.8778D0 
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4.94D+01 
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N ISUP HREF FIRST CARD 
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Table D-3. Input data for example given in Program 
Listing and Sample Output. 
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too severe, the user may modify the equations in subroutine F under a 
particular Key number to provide for any variation in power with angle of 
attack. Finally, the program structure is such that during one execution 
the same variable may not be specified in two different ways, since this 
would necessitate two different FUNCTION subprograms for the same calling 
statements. However, a variable may be incremented for several similar 
trajectories by simply passing a counter in the argument list of the 
FUNCTION call. For most light aircraft studies the aforementioned restrict- 
ions pose no problem to the user. 

The following discussion contains "trouble shooting" aids for the user 
who has a new problem for the program to solve. If the program message 
reads: 

(1) "INCONSISTENT INPUT DATA, EXECUTION CEASED FOR THAT KEY NUMBER" - 

Lift and weight were not equal or the program attempted to 
calculate a Cq less than Cq q at the initial point. 

(2) "VALUE OF MAXHLV WAS EXCEEDED, EXECUTION STOPPED" - The value of 

MAXHLV should be increased and/or use a smaller initial step 
size (TDEL) . 

(3) "ABSOLUTE VALUE OF FLIGHT PATH ANGLE BECAME GREATER THAN 1.0 

RADIAN" - If this error suddenly appears when the actual printed 
time history is small in magnitude but oscillatory, the corrective 
action if to decrease HMAX so as to prevent the rapid halving and 
doubling of the step size. 

(4) "ALTITUDE BECAME NEGATIVE" - This may occur when the actual tra- 

jectory closely approaches sea level. If the initial altitude 
is zero, then the start of integration may cause a slight oscil- 
lation which gives an erroneous negative altitude. In this case, 
increase h Q to perhaps 50 feet. 

For maximum integration efficiency, values for ACCHLV and ACCDBL may be 
adjusted based on the experience of previous executions to produce a 
minimum of halving and doubling the step size. 

The following listing and sample output complete this description of the 
Path Performance Program. 
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Sample Output 


POWER AVAILABLE. VS. VELOCITY 
REFERENCE ALTITUDE ■ 0.0 FEET 


KFT-LBS/SEC) 

VACFT/SECI 

0.0 


0.0 


0.29150D 

05 

0.273300 

02 

0.924 TOO 

09 

0.546700 

02 

0.699400 

05 

0. 820000 

02 

0.81620D 

05 

0.I0933D 

03 

0.874500 

05 

0.136670 

03 

0.909480 

05 

0.164000 

03 

0.944460 

05 

0.191330 

03 

0.958650 

05 

0.218670 

03 

0.967780 

05 

0.246000 

03 

0.973610 

05 

0.273330 

03 

0.979440 

05 

0.300670 

03 

0.994700 

05 

0.328000 

03 

0.99470D 

05 

0.35533D 

03 

0.994700 

05 

0.362660 

03 


HI • 

10000. 

VI - 

149.0 

GAN1 - 0.0 

ALPHA1 

- 5.878 

Ml - 

2650. 

PI • 

49.40 

XI - 

0.0 

C - 0.6000 

TOEL - 

0.5000 

TMAX * 

480.0 

6 - 32.20 

s » 

174.0 

RHO * 0.23600-02 

CLAO * 

0.3090 

CLA • 

4.608 

COO - 

0.26900-01 

E - 

0.9800 

AR - 7.378 

KEY - 

3 

IPRtNT 

• 10 

NAXHLV 

• 2 

ACCHLV • 0.10000-02 

ACCOBL * 0.20000-04 

HHAX - 

60.00 

VMIN - 

o 

e 

WENPTV 

- 2150. 

EX - 

2.000 







TINE 

ALTITUOE 


VELOCITY 

GANNA 

ALPHA 


CL 


CD WEIGHT 

POWER 

RANGE 

ININ! 

IFTI 


IFT/SEC 1 


IDES) 

I0EGI 




I POUNDS! 

(HP! 

(NILES! 

0.0 

10000. 



149.0 


0.0 

5.877 


0.7817 


0. 5380D-01 

2650. 

49.40 

0.0 

INTERVAL 

OOUBLEO AT 

T 


0.2500D 

01 

SECONDS 

NEW TOEL 


0.1000D 

01 

SECONOS ••*• 




•944 INTERVAL 

DOUBLED AT 

T 


0.75000 

01 

SECONDS 

NEW TOEL 


0.20000 

01 

SECONDS 44*4 




0.1250 

10000. 



149.0 


0.0 

5.877 


0.7817 


0.53800-01 

2650. 

49.40 

0.2116 

444* INTERVAL 

OOUBLED AT 

T 


0.1750D 

02 

SECONDS 

NEW TOEL 


0.40000 

01 

SECONOS 




4444 INTERVAL 

DOUBLED AT 

T 


0.37500 

02 

SECONOS 

NEW TOEL 


0.80000 

01 

SECONDS 4*** 




0.6250 

10000. 



149.0 


0.0 

5.877 


0.7817 


0. 53800-01 

2650. 

49.39 

1.058 

4646 INTERVAL 

OOUBLED AT 

T 


0.7750D 

02 

SECONOS 

NEW TOEL 


0.16000 

02 

SECONOS 4444 




4444 INTERVAL 

OOUBLED AT 

T 


0.1575D 

03 

SECONDS 

NEW TOEL 


0.32000 

02 

SECONOS 4444 




2.625 

10000. 



149.0 


0.0 

5.877 


0.7817 


0.53800-01 

2649. 

49.37 

4.444 

4444 INTERVAL 

DOUBLED AT 

T 

■ 

0.3175D 

03 

SECONOS 

NEW TOEL 

m 

0.60000 

02 

SECONOS 4444 




10.29 

10000. 



148.9 


0.0 

5.877 


0.7817 


0.5380D— 01 

2645. 

49.26 

17.42 

20.29 

10000. 



148.T 


0.0 

5.877 


0.7817 


0.53800-01 

2640. 

49.12 

34.33 

30.29 

10000. 



148.6 


0.0 

5.877 


0.7817 


0.53800-01 

2635. 

48.99 

51.22 

40.29 

10000. 



148.4 


0.0 

5.877 


0.7817 


0.53800-01 

2630. 

48.85 

68.09 

50.29 

10000. 



148.3 


0.0 

5.877 


0.7817 


0.53800-01 

2625. 

48.71 

.84.95 

60.29 

10000. 



148.2 


0.0 

5.877 


0.7817 


0.53800-01 

2620. 

48.58 

101.8 

70.29 

10000. 



148.0 


0.0 

5.877 


0.7817 


0.53800-01 

2616. 

48.44 

118.6 

80.29 

10000. 



147.9 


0.0 

5.677 


0.7817 


0.53800-01 

2611. 

48.31 

135.4 

90.29 

10000. 



147.8 


0.0 

5.877 


0.7817 


0.53800-01 

2606. 

48.18 

152.2 

100.3 

10000. 



147.6 


0.0 

5.877 


0.7817 


0.53800-01 

2601. 

48.04 

169.0 

110.3 

10000. 



147.5 


0.0 

5.877 


0.7817 


0.53800-01 

2596. 

47.91 

185.8 

120.3 

10000. 



147.3 


0.0 

5.877 


0.7817 


0.53800-01 

2592. 

47.78 

202.5 

130.3 

10000. 



147.2 


0.0 

5.877 


0.7817 


0.53800-01 

2587. 

47.65 

219.3 

140.3 

10000. 



147.1 


0.0 

5.877 


0.7817 


0.53800-01 

2582. 

47.51 

236.0 

150.3 

10000. 



146.9 


0.0 

5.877 


0.7817 


0.53800-01 

2577. 

47.38 

252.7 

160.3 

10000. 



146.8 


0.0 

5.877 


0.7817 


0.53800-01 

2573. 

47.25 

269.4 

170.3 

10000. 



146.7 


0.0 

5.877 


0.7817 


0.53800-01 

2568. 

47.12 

286.1 

180.3 

10000. 



146.5 


0.0 

5.877 


0.7817 


0.53800-01 

2563. 

46.99 

302.7 

190.3 

10000. 



146.4 


0.0 

5.877 


0.7817 


0.53800-01 

2558. 

46.67 

319.4 

200.3 

10000. 



146.3 


0.0 

5.877 


0.7817 


0. 53600-01 

2554. 

46.74 

336.0 

210.3 

10000. 



146.1 


0.0 

5.877 


0.7817 


0.5380D-01 

2549. 

46.61 

352.6 

220.3 

10000. 



146.0 


0.0 

5.877 


0.7817 


0. 53800-01 

2544. 

46.48 

369.2 

230.3 

10000. 



145.9 


0.0 

5.877 


0.7817 


0.53800-01 

2540. 

46.35 

385.8 

240.3 

10000. 



145.7 


0.0 

5.877 


0.7817 


0.53600-01 

2535. 

46.23 

402.4 

250.3 

10000. 



145.6 


0.0 

5.877 


0.7817 


0.53800-01 

2531. 

46.10 

418.9 

260.3 

10000. 



145.5 


0.0 

5.877 


0.7817 


0.53600-01 

2526. 

45.98 

435.5 

270.3 

10000. 



145.3 


0.0 

5.877 


0.7817 


0.53600-01 

2521. 

45.85 

452.0 

260.3 

10000. 



145.2 


0.0 

5.877 


0.7817 


0.5380D-01 

2517. 

45.73 

468.5 

290.3 

10000. 



145.1 


0.0 

5.677 


0.7817 


0.53600-01 

2512. 

45.60 

485.0 

300.3 

10000. 



144.9 


0.0 

5.677 


0.7817 


0.53800-01 

2506. 

45.48 

501.5 

310.3 

10000. 



144.8 


0.0 

5.877 


0.7817 


0.53800-01 

2503. 

45.35 

517.9 

320.3 

10000. 



144.7 


0.0 

5.87? 


0.7817 


0.53800-01 

2499. 

45.23 

534.4 

330.3 

10000. 



144.5 


0.0 

5.877 


0.7817 


0.53800-01 

2494. 

45.11 

550.8 

340.3 

10000. 



144.4 


0.0 

5.877 


0.7817 


0.53800-01 

2490. 

44.99 

567.2 

350.3 

10000. 



144.3 


0.0 

5.877 


0.7817 


0.53600-01 

2485. 

44.86 

583.6 

360.3 

10000. 



144.2 


0.0 

5.877 


0.7817 


0.53800-01 

2481. 

44.74 

600.0 

370.3 

10000. 



144.0 


0.0 

5.877 


0.7817 


0.53800-01 

2476. 

44.62 

616.4 

380.3 

10000. 



143.9 


0.0 

5.877 


0.7817 


0.53600-01 

2472. 

44.50 

632.8 

390.3 

10000. 



143.8 


0.0 

5.877 


0.7817 


0.53800-01 

2467. 

44.36 

649.1 

400.3 

10000. 



143.6 


0.0 

5.877 


0.7817 


0.5380D-01 

2463. 

44.26 

665.4 

410.3 

10000. 



143.5 


0.0 

5.877 


0.7817 


0. 53600-01 

2456. 

44.14 

681.7 

420.3 

10000. 



143.4 


0.0 

5.877 


0.7817 


0.53800-01 

2454. 

44.02 

698.0 

430.3 

10000. 



143.3 


0.0 

5.877 


0.7817 


0.53800-01 

2450. 

43.91 

714.3 

440.3 

10000. 



143.1 


0.0 

5.877 


0.7817 


0. 53800—01 

2445. 

43.79 

730.6 

450.3 

10000. 



143.0 


0.0 

5.877 


0.7817 


0.53800-01 

2441. 

43.67 

746.9 

460.3 

10000. 



142.9 


0.0 

5.877 


0.7817 


0.53800-01 

2436. 

43.55 

763.1 

470.3 

10000. 



142.7 


0.0 

5.877 


0.7817 


0.53800-01 

2432. 

43.44 

779.3 

480.3 

10000. 



142.6 


0.0 

5.877 


0.7817 


0.53800-01 

2428. 

43.32 

795.5 
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Hi * o.o 

VI - 

120.0 

GAM1 - 

1.500 

AIHHA1 - 7.420 

HI - 

2700. 

PI - 153.0 

XI - 

0.0 

C « 0 

» 

o 

o 

e 

TOEL - 0.5000 

TMAX - 

30.00 

G - 32.20 

S - 

174.0 

RHQ - 

0.23000-02 

Ct AO » 0.3090 

CLA - 

4.809 

CDO - 0.29870-01 

E ■ 

0.5930 

AH ■ 

7.379 

KEY - 12 

(MINT 

• 5 

MAXHLV - 2 

ACCHLV - 0.10000-02 

ACCD8L 

• 0.20000-04 

HMAX ■ 30.00 

VMIN - 

35.00 

WENPTV « 2150. 

EX - 

2.991 







TIME 

ALTITUDE 


VELOCITY 

GAMMA 

ALPHA 

CL 


CO 1 

HEIGHT 

POWER 

RANGE 

1 HI HI 

1 FT I 




1 FT/SEC 1 


( DEG) 

(DEG) 




(PCUNCSI 

(*PI 

(MILES! 

0.0 

0.0 




120.0 


1.500 

7.416 

0.9054 


0. 63920-01 

2700. 

153.5 

0.0 

*•44 INTERVAL 

00U8LE0 

AT 



0.25000 

01 

SECONOS 

NEH TOEL - 

0.10000 

01 

SECONDS 4*** 




0.41670-01 

8.223 




131.3 


1.500 

5.569 

0.7568 


0. *1490-01 

2700. 

157.4 

0.59490-01 

4444 INTERVAL 

DOUBLED 

AT 



0.75000 

01 

SECONOS 

NEH TOEL - 

0.20000 

01 

SECONDS 44** 




0.1250 

26.76 




151.7 


1.500 

3.208 

0.5670 


0.43200-01 

2700. 

162.4 

0.1937 

•••• INTERVAL 

DOUBLED 

AT 



0.17500 

02 

SECONDS 

NEH TDEL • 

0.40000 

01 

SECONOS 4*** 




0.2917 

70.76 




192.2 


1.500 

I. 049 

0.3934 


0. 34340-01 

2700. 

169.5 

0.5119 

444* INTERVAL 

DOUBLED 

AT 



0.37500 

02 

SECONOS 

NEH TDEL - 

0.80000 

01 

SECONDS 44*4 




0.6250 

175.7 




213.8 


1.500 

-0.2776 

0.2867 


0.3160D-01 

2699. 

172.9 

1.271 

44** INTERVAL 

DOUBLED 

AT 



0.77500 

02 

SECONDS 

NEW TDEL • 

0.16000 

02 

SECONOS 44** 




1.292 

410.3 




229.3 


1.500 

-C.7230 

0.2509 


0.3103D-01 

2698. 

172.5 

2.968 

4*4* INTERVAL 

00UBLED 

AT 



0.15750 

03 

SECONOS 

NEH TOEL - 

0.3000D 

02 

SECONDS 4**4 




2.625 

993.5 




231.0 


1.500 

-0.7283 

0.2504 


0. 3103D-01 

2696. 

169.6 

6.463 

5.125 

1799. 




229.8 


1.500 

-0.6170 

0.2594 


0. 31160-01 

2691. 

164.2 

13.01 

7.625 

2698. 




228.6 


1.500 

-0.4982 

0.2699 


0.31300-01 

2697. 

158.9 

19.52 

10.13 

3593. 




227.3 


1.500 

-C.3726 

0.2790 


0.3147D-01 

2683. 

153.7 

25.99 

12.63 

4483. 




225.9 


1.500 

-0.2394 

0.2897 


0.31660-01 

2680. 

148.7 

32.42 

15.13 

5367. 




224.4 


1.500 

-0.97870—01 

0.3011 


0.31890-01 

2676. 

143.8 

38.82 

17.63 

6245. 




222.9 


1.500 

0.52790-01 

0.3132 


0.3213D-01 

2672. 

139.0 

45.17 

20.13 

7116. 




221.1 


1.500 

0.2136 

0.3262 


0.3242D-01 

2669. 

134.3 

51.47 

22.63 

7981. 




219.4 


1.500 

0.3859 

0.3400 


0.32760-01 

2666. 

129.8 

57.73 

25.13 

8839. 




217.5 


1.500 

0.5710 

0.3549 


0.33I5D-01 

2663. 

125.4 

63.93 

27.63 

9689. 




215.4 


1.500 

0.7710 

0.3710 


0.33620-01 

2659. 

121.2 

70.08 

30.13 

0.10530 

i 05 


213.2 


1.500 

0.9881 

0.3885 


0. 34170-01 

2656. 

117.0 

76.16 


APPENDIX E - Lift-Drag Curve Fitting Program 

User Instructions 


Using a Least-Square-D i stance curve fit procedure (Ref. 26) this 
program yields a general drag polar of the form 

C D = k 1 + K 2 C 2 l + kjC^ 4 

where drag coefficient versus lift coefficient data is supplied. The program 
is written in FORTRAN IV and is designed to run in double precision on an 
IBM 370-165 computer with an average execution time of three to five seconds 
per curve fit investigated. Since this program is designed to be used in 
conjunction with the programs in Appendices C and D, the user has the option 
of four types of the general drag polar: 

2 k4 

(1) Cq = ki + k 2 CL + kjCL where all four coefficients kj, k 2 , k$, and 
k 4 are varied in the fitting process, 

2 k 4 

(2) Cq = Cq q + k 2 C|_ + kjCL where Cq q is specified by the user and k 2 > 
k-j, and k 4 are varied, 

k4 

(3) Cq = k^ + k 3 CL where kj, k 3 , and k 4 are varied. 

k 4 

(4) Cq = Cq q + k3CL where k^ and k 4 are varied and C Dq is specified 
by the user. 

For options (2) and (4) the first specified data point must be the zero-lift 
drag coefficient. Note that the functional form of the general polar 
prohibits the use of negative values of Cq. 

The program requires the specification of the following input data: 


(1) 


The number N of drag coefficient versus lift coefficient data 
points, the control parameter IKEY which specifies which of the 
four types of the general drag polar is to be used as the 
fitting functions. 


IKEY = 1 
IKEY = 2 
IKEY = 3 
IKEY = 4 


Cq = k*| + 

C D = C D n + k 2 C L + k 3^l. 4 


k 4 


C D = k 1 + k 3 C L 
C D = C D q + k 3 c L 4 » 


(2) The N data points of drag coefficient Cq versus lift coefficient 
C(_, one data point per card with Cq specified first in ascending 
order based on the magnitude of Cq. 
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Statements (1) and (2) above represent a complete set of data for one fitting 
process. Table 1 gives the format specification of the input data. Upon 
completion of the calculations with a complete set of data the program returns 
to the statement where N and IKEY are read. This gives the user the option 
of fitting another data set if desired. If N is positive when N and IKEY 
are read, a curve fit is performed on the new set of data. If N is zero 
the program terminates. Thus, the final data in the program must specify 
N as zero. 



Table E-1. Data format for curve fit program. 
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1 19 


DYOAL ( 3* I) - XX**ALI4) 

OVOAL(4 t l> - ALI3I*XX**ALI 4l*OLOC(XX) 

Y(I1 * AU1) ♦•AL l 2 l*XX**2 ♦ ALOI *XX»#AL < 4 » 



THIS SUBROUTINE SOLVES THE LINEAR SIMULTANEOUS EQUATION 


Sample Output 


DATA TO BE 


CD 

0.470000-01 
0.50000D-01 
0* 540000-01 
0.970000-01 
0. 600000-01 
0.700000-01 
0.800000-01 
0.100000 00 
0.120000 00 
0.140000 00 
0.162000 00 
0*162000 00 
Q. 20100D 00 
0.226000 00 
0.240000 00 
0.260000 00 
0.280000 00 
0.300000 00 
0.32000D 00 


FITTEO 

Ci 

0.0 

0.229000 00 
0.400000 00 
0.500000 00 
0 .570000 00 
0.730000 00 
0.830000 00 
0.970000 00 
C. 106500 01 
0.113500 01 
0.119000 01 
0.123000 01 
0.126000 01 
0.129000 01 
0.130500 01 
C.132000 01 
0.133500 01 
0.135000 01 
0.136000 01 


RNS - 
ALIM 

0.23450D-01 

—0.23614 

0.26216 

3.0000 

RNS * 
ALI1I 

0.176440-01 

—0.385730-01 

0.985200-01 

3.4610 

RNS - 
ALU! 

0.130550-01 

—0.514780-01 

0.11010 

3.6160 

RNS - 
ALIII 

0.205230-01 
• 0.655460-01 

0.214400-01 

5.1761 

RNS * 
Atm 

0.850840-02 
• 0.8787*0-02 

0. 449220-01 

5.3366 

RNS - 
Aim 

0.702210-02 
- 0.407890-01 

0.203320-01 

6.9554 

RNS - 

Aim 

0.519260-02 
« 0.307680-01 

0.222460-01 

7.1130 

RNS • 
ALII) 

0.354140-02 
» 0.420450-01 

0.134500-01 

8.4475 

RMS • 
Aim 

0.43574D-02 
• 0.463450-01 

0.876130-02 

9.5789 

RMS • 

Aim 

0.244940-02 
- 0.462370-01 

0.903850-02 

9.6964 

RNS « 
ALIII 

0.226590-02 
• 0.466310-01 

0.740080-0. 

10.303 

RNS - 
ALIM 

0.225110-02 
• 0.494730-01 

0.685760-02 

10.535 

RNS - 
ALIII 

0.225110-02 
- 0.494600-01 

0. 685400-02 

10.537 

RNS « 

ALIII 

0*225110-02 
• 0.494800-01 

0.665400-02 

10.537 

RNS * 
ALIII 

0.225110-02 
■ 0.494800-01 

0. 685400-02 

10.537 

RNS - 
ALU 1 

0.22511D-02 
* 0.494800-01 

0.68540 D-02 

10.537 

RNS - 
ALIII 

0.225110-02 
• 0.494600-01 

0.665400-02 

10.537 

RMS - 
ALII) 

0.225110-02 
• 0.494800-01 

0.665400-02 

10.537 

RMS - 

Aim 

0.225110-02 
■ 0.494600-01 

0.665400-02 

10.537 

RMS - 

Aim 

0.225110-02 
• 0.494600-01 

0.685400-02 

10.537 


********** AFTER 100 ITERATIONS 11.0 - RMS/RNSOID I • THE RELATIVE ERROR BETWEEN 

TWO SUCCESSIVE VALUES OF RNS, STILL EXCEEOS THE SPECIFIED ERROR PARAMETER - 


0.5000-02 


CD - 0. 470000-01 ♦ 


RNS « 0.225110-02 

0.494800-01*CI**2 ♦ 0.665400-02*0.** 0.105370 02 


APPENDIX F - Power Estimation 


A step by step procedure for estimating aircraft power available is well 
outlined in an appendix in Reference 7. While this procedure is very 
useful in predicting the power available in general, the maximum power 
available or the maximum power available for continuous operation may be 
found by using a simplified procedure. 

The type of propeller charts given in Figure F-1 has been the standard 
NACA design chart since 1929. The charts a I ready exist for many propellers; 
those having R.A.F. 6 and Clark Y airfoil sections can be found in References 
7 and 8. 


In sizing the propel ler for a new design, the first step after choosing 
a blade section is to calculate Cs from the equation: 

(0.638)(V n )(o) 1/5 

C_ = — 1 - — (F-1 

(BHP) 1/S (N) 2/5 

where 

Vq = design speed in mi les per hour, 
a = p/Pg = density ratio, 

BHP = design brake horsepower, 

N = design engine revolutions per minute. 


Using this value of C s project upward on Figure F-1 to the broken line 
of maximum efficiency for C s . This point determines the design blade angle 
and a horizontal projection to the V/nD scale gives the design value of V/nD 
where V = the velocity in feet per second, n = engine revolutions per 
second, and D = propel ler diameter in feet. The blade diameter can then 
be found by: 


( Vp) (88) 
(N) (V/nD) 


(in feet) . 


(F-2) 


The design efficiency is obtained by projecting upward from the design C s 
to the efficiency curve for the particular design blade angle. 

If the aircraft already exists the process of estimating the power 
available as a function of speed is somewhat different. One already knows 
the propeller diameter, the propeller section, the engine BHP at rated N, 
the maximum value of N for continuous operation, the type of propeller 
(fixed p i tch or constant speed), and the blade angle at 75 percent out the 
propeller radius. The maximum power available as a function of velocity 
is then calculated according to the procedures described below. 


Fixed Pitch Propellers 

The maximum horsepower for continuous operation at sea level is taken 
to be: 
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( F-3) 


(HP) = ( BHP ) ( N /N )(n) . 
max c r 

where 

BHP = engine brake horsepower at N r , 

N c = maximum engine revolutions per minute 
for continuous operation, 

N r = engine rated revolutions per minute, 
r) = propeller efficiency. 

The propeller efficiency for a given blade angle at 75 percent of the pro- 
peller radius is a function only of (V/nD) or V since n(n = N c /60) and D 
are specified. Choosing values of (V/nD) and using the efficiency curves 
for the particular blade section, the value of ( p P) m and V can be calculated 
for each value of (V/nD) chosen. One would thus be able to complete a table 
similar to the one below. 


(V/nD) 

(HP) v 
max 

from Equation (F-3) 

V 

V = (V/nD) ( nD) 

n 

(from section chart) 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


Table F-1. Sample tabulation of maximum power and velocity. 


Engine efficiency, if known, can be used to multiply the right hand 
side of Equation (F-3) so as to give a more accurate value of (HP) max . 

The power available at any altitude is obtained by multiplying the sea 
level power required by the factor (a - 0.165)/0.835 where a = p/p Q = the 
ratio of the density of air at altitude to the density of air at sea level 
(formula obtained from Reference 3). 


Constant Speed Propel lers 

The maximum horsepower for continuous operation at sea level for a 
constant speed propeller is also given by Equation (F-3). Choosing values 
of (V/nD) efficiency values can be obtained from the envelope of the 
efficiency curves since the best efficiency at a particular (V/nD) may be 
achieved by varying the propeller blade or pitch angle. Since n and D are 
constant for the maximum power case, one can find n for any value of V. 
(HP) ma x f° r an y velocity is then readily determined . The results may be 
tabulated as in Table F-1. 
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APPENDIX G - A Predictor Corrector Method for Numerical Integration 

of Relaxation Differential Equations 


By Neill S. Smith 


Relaxation differential equations are character i zed by a strong depend- 
ence of the derivative of the dependent variable on the difference between 
its own value and a slowly varying function. Conventional Runge-Kutta and 
predictor-corrector methods are unable to handle these types of equations. 
Generally, these conventional methods develop strong oscillations when 
applied to relaxation differential equations. 

Trenor (Ref. 31) derived a modification of a fourth-order Runge-Kutta 
method for use with relaxation differential equations. His method is 
particularly useful because it becomes identical to the conventional Runge- 
Kutta in regions where the derivative is not strongly dependent on the value 
of the dependent variable. Thus his method can be applied to a differential 
equation which is not of relaxation form over the entire integration range. 

Controlling the integration step size in order to maintain a specified 
accuracy is difficult with Runge-Kutta methods. The usual practice is to 
integrate from the point x c to the point x Q + h using the step size h and 
then integrate again from x Q to x Q + h using the step size h/2. The dif- 
ference between the two results obtained at x Q + h is used to determine 
whether the step size should be halved, doubled, or remain the same. This 
procedure requires eleven evaluations of the derivative function to integrate 
forward one step. 

On the other hand, controlling the step size to maintain a specified 
accuracy is fairly simple with predictor-corrector methods. The difference 
between the corrected and predicted values or the difference between two 
successive interated values of the corrector can be used to control the 
step size. Thus the pred ictor-corrector method requires only two or three 
evaluations of the derivative function per integration step. 

Obviously there results a great savings in computational time with 
predictor-corrector methods; therefore, a pred ictor-corrector method that 
could handle relaxation differential equations would be extremely useful. 

Such a predictor-corrector method is derived in this Appendix by applying 
Trenor's modification to the corrector of the conventional Adams-Bashforth 
predictor-corrector method. In regions where the differential equation is 
not of relaxation form, the modified corrector becomes identical with the 
conventional Adams-Bashforth corrector. 


The procedure for integrating the first-order differential equation 


t 

y 


- f(x ^ 


(G-l ) 


from the point x n to x n +i = x n + h by the conventional Adams-Bashforth 
method is given below: 
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(1) An estimate of the value of y at the point x n+ ^ denoted by m n+ ^ 
is obtained with the predictor equation 

m n+1 - Vn + 24 [55 Vn ' 37 Yn-1 + '5 Vn- 2 - 9 Yn-j] • <G - 2 > 

(2) An estimate of the derivative of y at x n+ ^, denoted by m n +i, is 
obtained by evaluating Equation (G-1) at the point (x n+1 , m n+ j). 

m n+1 = f(x n +1 » m n+1 > (6-3) 

* 

(3) Using m n +i (the estimated derivative of y at x n+ ^ ) a final 
corrected value of y at x n+ ^ is obtained with the corrector 
equation 

¥ n+1 = Yn + ^9 m n+1 + 19 y n - 5 y^ + y n _ 2 ] . (G-4) 

Although generally very satisfactory, this method fails when Equation (G-1) 
is of relaxation form, i.e. , it can be written approximately as 

^ = - P(y - y) (6-5) 

% 

where P is a large positive number, and y is a slowly varying function of x. 

A modification of the corrector Equation (G-4) that enables the Adams- 
Bashforth method to handle the above situation is derived below. 


Following Trenor's procedure, it is assumed that Equation (G-1) can be 
approximated by 

= f (x,y) = - P(y - y ) + A + B(x - x ) + C(x - x ) 2 + D(x - x„) 3 (G-6) 

dx n n n n 

over the interval from x n _3 = x n - 3h to x n+ ^ = x n + h. Equation (G-6) can 

be integrated to give 


y n +1 = y n + H [AF 1 + hBF 2 + 2h 2 CF 3 + 6h 3 DF 4 ] 
where the functions F n are simple exponential functions of Ph 

F q = exp [-Ph] 

F = F n-1 ~ (n-D! _ “ (-Ph) k 

n (-Ph) k=0 (n+k)! 


(G-7) 


(G-8) 


The five constants A, B, C, D, and P are evaluated by requiring Equation (G-6) 
be satisfied at the five points x n +j, i = 1 ,0,-1 ,-2,-3. The resulting 
expressions for these five constants are 
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= y 

Hzl _ z + fn + !n+i 
6 n-1 2 3 



2h 2 C 


z n-1 ‘ 


2Z„+Z 


n+1 


(G-9) 


bh 3 u - - Z n _ 2 + 3 Z n _ 1 - 3 Z n + Z n+1 

P = - Vn+1 ~ 4 ^n + 6 Vn-1 ~ 4 Vn-2 + Vn-3 

Vn+1 - 4 Vn + 6 Vn-1 " 4 y n -2 + / n -3 

t i » 

where Z n+ j = y n +; + Py n+ j. Since y n+ ^ and y n +i are used to determine the 
five constants. Equation (G-7) represents a corrector equation that is to 
be used in place of the conventional corrector given by Equation (G-4) . When 
P = 0 (corresponding to no relaxation-type dependence of y on y), then = 1, 
F2 = 1/2, F3 = 1/6 and F4 = 1/24, and Equation (G-7) becomes identical to the 
conventional corrector given by Equation (G-4). Since P is continually 
calculated as the integration proceeds, the predictor-corrector method obtained 
by replacing the corrector given by Equation (G-4) with the modified corrector 
given by Equation (G-7) will automatically handle any relaxation dependence 
that appears in the differential equation. In regions where the differential 
equation has no relaxation dependence (P = 0), the method automatically 
becomes identical to the original predictor-corrector method. 
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APPENDIX H - A Discussion on Specific Fuel Consumption 


Throughout this work it has been assumed that the specific fuel 
consumption, c, of a piston engine is constant. This of course is a rather 
gross approximation and it is the purpose of the present section to 
determine how serious the error is, how one may incorporate a more exact 
model in the calculations if desired, and how the use of a more exact 
model will alter the conclusions reached previously. 

A piston engine is a very complex machine; thus, unless one uses 
experimental data giving power and fuel flow rate as functions of manifold 
pressure and shaft speed in a table look-up form, it is necessary to make 
some approximations to these characteristics in order to obtain managable 
functional forms. The experimental characteristics shown in the figures 
below indicate the magnitude of the problem. 


BHP 



Figure H-1 . 


Typical variation of BHP with RPM for 
various constant throttle settings. 
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settings. 
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H 1 


1-3. Typical variation of SFC with BHP for 
various constant RPM settings. 
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The constant speed contours in Figure H-3 are given for each 100 RPM 
above fast idle. 

In terms of the throttle position, it, and RPM, N, a fairly good 
representation of the power output is 



where the value of Cq and C 2 depends upon the particular engine involved. 

During most of the flight regime one could normally attempt to operate 
at minimum fuel consumption. On an aircraft equipped with a constant speed 
propel ler this means one will generally move down a constant N contour 
(Figure H-3) until full throttle is reached and then along a constant 
manifold pressure contour. This type of operation can usually be represented 
fairly accurately by a function of the type 

c = + A 2 P 1 2 3 . 

P z 

Obviously, more complex functions can be used to obtain more accurate 
descriptions of experimental results. 

For the moment suppose the foregoing expression for c represents the 
physical situation adequately. How then does one incorporate it in the 
performance analysis, and how are the previous conclusions one obtained 
on the basis of the simple, linear f ue I -f I ow-power function altered? For 
the more general fuel consumption function, 

W = A 1 + A 2 P 4 . 

Mechanically, such a function presents no substantial difficulties in 
the computation procedure, although the program instructions would have 
to be changed to accomodate it. Generally, it is to be expected that 
such a fue I -f low-power relation does not alter the speed for maximum 
endurance and reduces the range achievable with high power (>75$) 
trajectories . 

If one wished to perform the computation precisely, the following 
procedure is suggested: 

1. Obtain from the engine manufacturer the most reliable data he 
has on the engine in question and plot the following curves: c versus P 
with N as parameter, c versus P with manifold pressure as parameter. The 
plots should include the effects of a particular propeller’s efficiency 
so that power is really thrust horsepower. 

2. Employ a least squares polynomial fit to each family of curves. 

3. Then when one specifies a manifold pressure and an engine speed, 

he solves the two equations simultaneously to obtain c and P. 
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4. The path performance computation is changed as follows: 

(a) Either manifold pressure or engine speed or both.must be 

specified. If both are specified c and P are both found and W can be 
computed. If only one is specified the power required which comes from 

a solution of the system of equations in Appendix B must be used in 
conjunction with the two curve fits to find c and engine speed or 
manifold pressure, whichever was previously unknown. 

(b) The value of c corresponding to the required value of P 
is then inserted in the path performance equations and the computation 
continues for another small increment in time. 

The magnitude of the programming task is quite evident. Unfortunately, 
the pressures of time did not permit its completion during the preparation 
of the present work. 
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